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THE TRANSMISSION OF EXCITATION FROM 
THE MEMBRANE TO ACTOMYOSIN 


Z. Bay, M. C. Goopatt, AND A. SzENT-GyORGYI 


INSTITUTE FOR MuscLE RESEARCH 
MARINE BIOLOGICAL LABORATORY 
Woops Hore, MASSACHUSETTS 


The impressed field, ‘‘Window Field’? (WF), due to a half-wave action potential on a 
muscle fiber, has been calculated on the basis of potential theory. It has been shown that in 
spite of the small intensity of the field, its integrated action can transfer the energy needed to 
induce contraction from the membrane to the interior of the fiber. The energy of polarization 
has been found to be sufficient to exceed the energy of thermal agitation on that length of 
fiber, which can be identified as the length of a sarcomere. The changes of ion concentration, 
caused by the WF, if calculated on the assumption of the semipermeability of the Z mem- 
branes, was found to be equal to the changes necessary to induce contraction of actomyosin 
in vitro. 


The contractile matter in cross-striated muscle is a complex protein, 
actomyosin (Szent-Gyérgyi, 1951). At rest the two component proteins, 
actin and myosin, are kept apart by coulombic repulsions. In vitro a 
slight change in charge, brought about either by a decrease of 0.01 M in 
the salt concentration, of 0.0001 M in the ATP (adenoside triphosphate) 
concentration, or a pH change of 1.5 units, may suffice to bring the two 
proteins together to form the actomyosin complex, which in the presence 
of ATP goes over spontaneously into its contracted state. Im vivo an 
analogous disturbance is brought about by “excitation.” The problem is 
whether depolarization of the membrane occurring in excitation can gen- 
erate the disturbance of charges necessary for actomyosin formation. 
This problem has been dealt with in a previous note (Bay and Szent- 
Gyérgyi, 1951). The present paper gives a more detailed treatment of the 
problem. 

Since the muscle fiber is surrounded by a homogeneously polarized 
membrane, every point in its interior has the same electric potential and 
there is no electric field. Obviously, if part of the membrane is depolar- 
ized, this situation is altered and the change can easily be calculated by 


simple electrostatics. 
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In the upper part of Figure 1, the two parallel horizontal lines represent 
the longitudinal section of a fiber, the membrane of which is depolarized 
to the right of z = 0 (the depolarization being symbolized by the dotted 
line). Before depolarization occurred, the potential at any point in the 
interior of the fiber was 47m = constant; m, the density of the dipole 
moment of the membrane, being constant. The partial depolarization of 


the membrane acts as if a “window” had been cut in the membrane. The 
potential at any point will now be 


V =m (4x — fdQ) =m (4r—Q), (1) 


if Q is the solid angle of the window. Since the window is a circle, its 
solid angle can easily be calculated at any point along the axis: 


Substituting (2) into equation (1) we find 
z 
a Voa(1— eR) (3) 


Vo being the resting potential of the membrane. 
The intensity of the electric field along the axis is E, = —dV/dz and 


Vo fs 


Ee 7R RT By (4) 
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The situation is represented in the lower part of Figure 1, which shows 
the cross-section of the fiber with its dipole layer maintained on the left 
and depolarized on the right side of s = 0. The potential is shown by the 
sloping line (circles). It changes inside the fiber gradually from V» to 0.* 
Its gradient is shown by the other full line (triangles). The maximal in- 
tensity of the field is 
Vo 


E, max — OR . 


This intensity would be obtained by setting the potential difference Vo 
on the plates of a condenser separated by the distance 2R, which is equal 
to the diameter of the fiber. This field will henceforth be called the 
“Window Field,” (WF). Taking for Vy and 2R the plausible values 0.1 
volt and 0.05 mm. we find 


E, max — 20 volt/ cm. 


The above presupposes that the depolarization of the membrane takes 
place with a sharp border. If this is not the case, the extension of the WF 
along the axis will be greater and its intensity smaller. However, it is im- 
portant to note that for the WF, propagated with the uniform velocity 2, 
the time integral of Z, (during the half-wave of depolarization) is inde- 
pendent of the sharpness of the border between the polarized and de- 
polarized parts of the membrane, i.e., is independent of the form of the 
propagated wave. To show this we take into account that in the case of 
such a propagation parallel to the Z-axis 


tin cpa 3 
oe Sore (>) 
and we have for the time integral 
a2 be Gare V Vio—V_o Vo 
= — — = — = -, 6 
SJ, Fatt nigh es ae v v 2 


Equation (6) also shows that the time integral of £, is identical for any 
point of the cross-section. Thus the time integral of the flux of the electric 
field through the cross-section is (V o/2)f, if f is the area of the cross-section, 
and it is independent of the shape of the cross-section of the fiber. 

The time integral of the component of the field vertical to the axis is 0 
for any point of the cross-section. 

* If the potential is taken to be 0 inside the polarized tube, the potential curve in Figure 1 


should be shifted downward by Vo. The gradient of its field will remain unchanged. In case of 
reversal of membrane potential during activity, Vo means the amplitude of the action potential. 
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To show this we now use another description of the WF, deriving it 
from a vector potential, A. That is, we apply the method used for the 
calculation of the field of magnetic dipole layers, surrounded by electric 
currents (Biot-Savart rule). In this case, the tangent of the border line 
between polarization and depolarization is perpendicular to the Z-axis 
and therefore ay (7) 
We consider one point in the plane of Figure 1 and adapt the coordinate 
system with the Y-axis in the plane of Figure 1, the X-axis perpendicular 
to it. 

Then we have 


and 


(% a=+ f eee, ae (8) 


The derivation of equation (8) shows that it is true for any form of the 
propagated wave and it is-also independent of the form of the cross- 
section of the fiber. 

The component of the WF, perpendicular to the axis, thus acts as an 
alternating current, depolarizing the membrane in front and repolarizing 
it behind the wave. This tends to make this latter self-propagating, if the 
membrane has the type of negative impedance which has been treated 
theoretically by G. Karreman (1951). 

In view of the large interionic fields (~ 10’ volt cm.) it can be expected 
that the quantity of actomyosin brought together depends not on the 
instantaneous intensity of the WF, which is quite negligible (~1 volt 
cm.—' for a wave of depolarization with un-sharp border), but on its inte- 
grated action which depends on Vo. In the following calculations we will 
show that such integrated action of the WF can transfer the energy needed 
to induce contraction from the membrane to the interior of the fiber. 
This seems to be of great importance, owing to the fact that the WF ap- 
pears and acts on the entire cross-section of the fiber without time delay. 

Before going into details we may consider the energy available in the 
membrane in relation to the possible elementary structures in the fiber. 
The membrane energy per unit length is 


Wm = RCV? , (9) 
where C is the capacity per unit area. This sets an upper limit for the 


energy available to change the equilibrium inside the fiber during the 
process of excitation. 
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The fiber as noted below consists of N filaments per unit area of cross- 
section, 1.e., the total number of filaments inside the fiber is R27. The 
stability of the system against thermal motion suggests that the energy 
W, needed to trigger the system is greater than R’xN&T, or it is 


W,=mnR?ax NRT , (10) 


where 7 > 1 and, say,  ~ 3. This does not imply that the total number 
of degrees of freedom in such a macromolecule is one or near one, but it 
seems plausible that one degree only is involved in bringing actin and 
myosin together and thus triggering contraction. We use m as a safety 
factor. 

Then assuming that the energy of the membrane will be transported 
by some reasonably efficient process to the filaments, one can ask for the 
length / of the fiber which has sufficient energy on the outside membrane 


to trigger the inside. This requirement gives us, using equations (9) 
and (10), 


<= : (11) 


One finds for muscle fibers C ~ 5 X 10-* farad cm.—*, Vo = 107 volt. 
With regard to the experimental evidence, the electronmicroscope reveals 
that the contractile matter consists of thin (~ 200 A wide) filaments 
running parallel to the axis, as shown first by C. E. Hall, M. A. Jakus, 
and F. O. Schmitt (1946). So V = 4 X 10" and we have 


l a 


It is striking that muscle, as seen in the microscope, actually has repeat 
units (sarcomeres of the length /, ~ 2.5 X 10-* cm.) of the same ratio 
1,/R ~ 10 as calculated here. This seems to show that the cross- 
structures must be considered in the mechanism of excitation and, in 
fact, may be the most significant part of it. 

The problem now is: How is the WF capable of this energy-transfer? To 
answer this question we calculate some integrated actions of currents of 
the WF, producing changes of concentration during the half-wave of de- 
polarization. 

In these calculations we have to use the longitudinal component (Es) 
of the WF, since it can be shown that the action of the transversal com- 
ponent is negligible. It was stated previously [eq. (8)] that the time inte- 
gral of the transversal component of the WF is 0 for the whole wave of 
depolarization. However, since the transversal component is an alternat- 
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ing current the value of its time integral is 0 for its two different phases, 
and one could ask whether either one of them could elicit contraction? 
That the time integral for one phase is negligible can be shown by the 
following calculations. 

The straightforward way of calculation would be to evaluate (8) be- 
tween the limits — © and the time at which Z, changes sign and at which 
the vector potential reaches its maximal value. This, however, would 
meet mathematical difficulties and for this reason we use here an approxi- 
mative method starting from the equation of divergence which is 0 inside 
the fiber. 

Let us denote the cylindrical coordinates of a point inside the fiber 
by p and z (p the distance of the point from the axis); then the equations 
of divergence are 


eV ,10V , eV _ 
Vii= ea, aS (13) 
and 
O02 Vien One 02V 
Pat Op Pras (14) 
or 
fe) OV 0? V 


For the present approximation we take 0?V/dz? independent of p (the 
longer the wave of excitation the better is this approximation) and inte- 
grate (15), obtaining 


Paes aaeeadal (16) 
or 

OV p 0?V 

5 Guaicene Oh etal (17) 


—0V/dp is the radial (transversal) component of the WF (the former £,) 
and its time integral is 


hav, pt rt avV 1/aV 
— f'—ai=2- dz=S-( ise (18) 


where fy is the time at which dV /dp changes its sign and at which dV /dz 
reaches its maximum value for the coordinate sz. 

Equation (18) shows explicitly that contrary to the time integral of the 
longitudinal component, the time integral of the transversal component 
depends on the sharpness of the border between polarization and de- 
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polarization because 0V/0z diminishes if the wave length increases. Let 
us use as the order of magnitude for (0V/dz) max the quantity V /1,, where 


14 is the length of the wave of excitation, and we obtain for the ratio of 
(18) and (6) 


1 p 
ee (19) 


Since in the actual case p < ~, mm. and/, ~ 1 mm., the ratio of the two 
time integrals is ~ 10-*. Thus in the subsequent treatment we can neglect 
the transversal component of the WF in our calculations, and we can 
base these calculations exclusively on the longitudinal component. 

On the subsequent pages we will consider electric currents inside the 
fiber due to the WF, i.e., the following treatment will be considered as a 
perturbation of the former equations. The problem will be whether these 
currents can cause changes in concentration of ions which can induce 
contraction. 

In a homogeneous and isotropic medium the WF cannot produce space 
charges or differences in concentration, because the WF is divergence- 
free inside the fiber, that is, 


V:-E=V:VV=0. (20) 


Introducing the current density 7 we have 
AS. (21) 


where ga is the conductivity and 


V-1=V- (cE) =oV-E+EVo. (22) 


Equation (22) shows that V - 7 can be different from zero if o is inhomo- 
geneous or anisotropic. 

Up to this point we have used only well-established facts of biology, 
such as the existence of a cylindrical membrane surrounding the fiber, 
polarized in the resting state and depolarized in excitation; therefore the 
above results are straightforward consequences of the potential theory. 
In the following we will give an example of a possible application of the 
theory outlined above. In order to do this we introduce an assumption 
which is connected with the well-known fact that, as mentioned above, 
the muscle fiber shows along its axis inhomogeneities (cross-striations) 
repeated with the periodicity /,, as, for example, the Z-membranes. We 
will suppose that these are less permeable for positive ions than for 
negative ones. These Z-membranes are not the only such inhomogeneities 
of the same periodicity. There are also M-membranes, I-discs, and Q-discs, 


8 Z. BAY, M. C. GOODALL, AND A. SZENT-GYORGYI 


the latter seeming to have a greater density and accumulating certain 
constituents of muscle, like glycogen and ash, as shown by A. B. Mac- 
Callum (1905). 

Due to these inhomogeneities, it is reasonable to assume that the con- 
ductivity is a function of the axial distance z, as shown in F igure 2. The 
distance J, is the thickness of a cross-membrane (e.g., the Z-membrane) 
which we suppose to be less permeable to positive ions (for the sake of 
simplicity we take here o+ ~ 0 inside the interval 1,*) and permeable to 
negative ions (i.e., ¢~ is the same inside and outside). 


FIGURE 2 


Applying the impressed field EZ, (henceforth simply £) the intensities 
of the electric field will be different inside (£;) and outside (£,) the cross- 
membrane. The total electric current is divergence-free, i.e., 


Exo =F, (ot+oa7). (23) 


The separate currents of the positive and negative ions, however, are not 
divergence-free. We obtain on one side of the membrane accumulations 
of positive and negative charges of equal amount, which are given by 


E,a0+ = (EE, —E2) o— (24) 


for unit area of cross-section and unit time. This means that an increase 
of the total concentration of the solution occurs during the flow of cur- 
rents. The same holds true with change of signs for the other side of the 
membrane, where the concentration decreases during the flow of currents. 

These changes of concentration, dependent on the integrated action 
of the WF, represent a useful transfer of energy from the surrounding 
membrane to the inside of the fiber during the half-wave of depolarization. 


* But it is not essential for the following calculations to have «+ = 0 exactly. 


TRANSMISSION OF EXCITATION 9 


The calculation of £, and Fy, as functions of time and dependent on £, 
involves the calculation of the different potential differences inside and 
outside the cross-membrane and the counter electromotive forces of ihe 
diffusion-potentials. We restrict ourselves to some simplified calculations 
which give the orders of magnitude. 

The ratio E,/, is independent of the impressed field and of time, and 
is determined by (23) to be 


a — 
Ey _¢ oe aE samhats (25) 
F2 o a 


using 


At the beginning of the process the diffusion potentials are negligible 
and we have the uniform fields E; and Fp, for which 


Fil, + Fol, =E (1,+1,) . (26) 
From equations (25) and (26), using 8 = /./1,, we obtain 


pite £46 
a (Gn eee 


a 


| 


(27) 


and 
US cats 


fag 


f,=E (28) 


a 


From the electronmicroscopic picture we have 


1 


B~55 


for the Z-membranes, so we can use the approximations 
i+a 1 (29) 
a 


Ky, =E 
ie 
a 


and 
E,~E Ay ‘ (30) 
1+= 
a 
We can estimate a from the measured average conductivity o of the 
fiber, writing 
Ec=E,(o++<0-). fea) 
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From equations (30) and (31) 


c= sed Rialglble Ss, 

fae 1+— 

a a 

so that 

a= z = — B (32) 

ot {Dex ? 

=—j|] —-—=—] 

o Dame: 


where 20+ can be taken as the conductivity of 0.15 M KCI. B. Katz (1948) 
found experimentally 20+/o = 3, so that 
os —~ 33 
a 75: (33) 
During the build-up process of the changes of concentration the ionic 
charges appear at the surface of the cross-membrane with the speed given 
by (24) and flow away by diffusion with the speed 
RI 2ut ae aC 2a sie 


Soe a SS ES ——ot 34 
c ut+tu dz 1l+a cdz ’ Sa. 


where u+ and wu are the mobilities of the + ions, V = kT/e, and c is the 
concentration. The maximum change in concentration occurs when the 
two speeds given by (24) and (34) are equal, i.e., 

da_ide 


eee 
= ia da 


ae (35) 


We use the notation ¢o for the concentration in the bulk of the solution, 
dc for the departure from ¢o at the surface of the cross-membrane, Jp for 
the distance which represents the region in Figure 2 where the concentra- 
tion change is appreciable, so that 


da det 
1l+a Co lp’ 


Emus * V (36) 
where we took E, ~ E and used the maximum intensity of the WF as 
that at which the maximum change of concentration occurs. 

We obtain another equation for 5c/cy and Ip, if we calculate the total 
charge passing per unit cross-section during the half-wave of depolariza- 
tion, determined by the time integral (6) of the WF. We have 


TRANSMISSION OF EXCITATION 11 


and 
p= ut, (37) 


We use the mobility w+ = 6 x 10-4 cm./sec./volt/em. of K+ and 
Emax ~5 volt/cm. given by the records of action potentials, V = 100 
cm. sec.—;so that 


6¢ 1 
he Spat —1 
= 10 (38) 
and 
lp ~2X10->cm.=2000A. (39) 


The value of 6c/co given in (38) is of the same order as the value ob- 
tained from experiments im vitro, being sufficient to bring actin and myo- 
sin together. Thus the WF can act on the length J, ~ 2000 A during the 
half-wave of excitation. 

To show that the results obtained in (38) and (39) are reasonable, we 
note that the counter diffusion potential per Z-membrane is of the same 
order as the impressed action potential on the length /,. We use Nernst’s 
formula: 


Se ee (40) 
i 


which we can easily apply outside of the cross-membrane and find 


- d¢cl1—a 
6V=2V Tau are (41) 
for the sum of the diffusion potential differences on the two sides of the 
membrane. Inside, the diffusion potentials are dependent on the value of 
a+ which we took as negligibly small compared to o+ outside, but which 
can be of the same order as o-, giving negligible diffusion potential dif- 
ferences in this case. 

Calculating only the value of (41) we obtain 


Oe? X10 volt (42) 
At the end of the process of the concentration polarization, the in- 


tensity of currents diminishes, so that the counter diffusion potential 
represents the only potential drop to be overcome by the WF, which gives 


El, ~ 10-3 volt (43) 


for one repeat unit. 
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The energy transferred by the WF during this process is also of the right 
value. To show this we calculated the isothermic work W done by the 
electric forces against the forces of diffusion while changing the concentra- 
tion on the two sides of the cross-membrane. The energy per unit cross- 
section is 

wer f Oe dx 6 Pcl, klic ae 1 erx~102RT, (44) 
Zz 0 


C Got os 


where cy = 0.15 M or cp = 107° ions/cm.’ and this is of the same order 
as wNRT, given earlier, showing that an efficient transfer of the membrane 
energy can occur through the WF. 

It is evident that the depolarization of the cell-membrane will revert 
the changes in ionic concentrations which the WF has brought about. 
However, this will happen only after the WF has done its work and the 
changes in ionic concentration have triggered off the contraction of the 
actomyosin system. Once actin and myosin have got together (in presence 
of ATP) the contraction-cycle is started up and will have to run to com- 
pletion. 

According to the theory outlined in this paper, it is the longitudinal 
field which brings the muscle to contraction. If this is true one has to 
expect to be able to bring the muscle to contraction by a longitudinal 
electric field of the calculated intensity of the WF, a few volts per cm. 
That muscle can be brought to contraction by an external longitudinal 
field was shown lately by W. A. H. Rushton (1930) who also reviewed 
the literature of this problem. He even showed that an external cross- 
field has no such effect. This experiment gives no information about the 
possible action of an internal longitudinal field. Rushton explained his 
results by supposing that the longitudinal external field brought the 
muscle to contraction by depolarizing the membrane at the cathode end 
of the muscle where the current leaves the fiber through the membrane 
in cross-direction. The local depolarization thus produced sets up a 
propagating wave which causes contraction of the muscle by running 
along its whole length. The length of the directly depolarized part at the 
end is of the order of 1 mm. 

In order to obtain evidence of the contracting action of an internal 
longitudinal electric field one has to place the muscle into a longitudinal 
field but at the same time abolish the propagation in the membrane and 
show that on application of the longitudinal field the muscle contracts in 


its whole length. As is generally known, the propagation in the membrane 
can be abolished by novocaine. 
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Our experiments were set up in the following way: similar to the ar- 
rangements of Rushton, the muscle was placed into a square trough filled 
with Ringer’s solution with two sheets of platinum at two opposite sides 
serving as electrodes. The sartorius muscle of Rana Pipiens was provided 
at both ends with ligatures at the tendon and suspended in Ringer’s 
solution under slight tension longitudinally to the field. The muscle was 
curarized and then 0.1% novocaine added gradually. It took 10-15 
minutes for the drug to abolish conduction. 

Square wave pulses were generated by a common type electronic pulse 
generator (Grass Instrument Co., Quincy, Mass.). In order to match the 
impedance of the trough (about 200 ohms) a cathode follower was used 
which was able to give output currents as high as 500 millamps. 

On applying the pulses the muscle started contracting when the in- 
tensity of the field reached the calculated intensity of the WF. The muscle 
contracted in its whole length and the intensity of contraction increased 
in the beginning with the intensity of the field and its duration. The 
rheobase (double of the minimum of the field intensity which produced 
reaction at long duration, >20 msec.) was around 1 volt per cm. 

It seems not impossible that the WF is also involved in more primitive 
forms of motion, as, for example, the propulsion of pseudopodia in pro- 
tozoa. If the membrane of a protozoa is depolarized locally at a point, the 
WF generated at this point must tend to eject negatively charged protein 
particles through the “window” from the interior, the work being done 
by the polarized part of the membrane which decreases this potential. 
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Some general properties of the solution of the diffusion equation are deduced for the steady- 
state, spherically symmetric system. On the basis of these developments some results of 
N. Rashevsky (Bull. Math. Biophysics, 11, 15, 1949) are discussed and the results of a previous 
investigation (Hearon, Bull. Math. Biophysics, 12, 135, 1950b) are extended to more general 
conditions. In particular these extensions apply to the flow of a solute against its concentration 
gradient, the nonzero gradient of an inert metabolite, and the accumulation or exclusion of an 
inert metabolite in a metabolic system. 


Introduction. lf there is considered a system in which the amount of 
a solute changes with time as a result of transport and production or con- 
sumption (due to chemical reaction) the equation of continuity for the 
kth chemical species is 


= aval 
ay ho Se, (1) 


where Te Q,, and C;, are the current (or flux), the rate of production, 
per unit volume, and the concentration of the kth substance, respectively. 

Nonlinearities may enter into (1) in two ways: The diffusion co- 
efficients, which appear in the J x, may depend upon the concentrations 
and the Q; may be nonlinear functions of the concentrations. The dif- 
fusion coefficients will, if the system is structurally nonuniform, also 
depend explicitly upon the space coordinates and the same may be true 
of the Q,. Indeed, the Q; may depend upon the space coordinates alone, 
as, for example, when they are functions of one or more catalysts, spatially 
distributed in the system, and every catalyst is saturated. A structural 
nonuniformity of the system which results in the explicit dependence of 
the diffusion coefficients upon the space coordinates does not necessarily 
imply that the Q; possess this property. 

* A portion of this work was performed while the author was a research participant, Oak 
Ridge Institute of Nuclear Studies, assigned to the Mathematics Panel, Oak Ridge National 
Laboratory. 
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In general (1) cannot be solved in the nonlinear case even for the 
steady-state condition. Therefore, methods which allow some informa- 
tion about the nature of the solution to be deduced from general proper- 
ties of the J ., the Q;, and the geometry of the system are of interest. 
Rashevsky (Joc. cit.) has discussed a steady-state case wherein it was as- 
sumed that the diffusion coefficient was constant and that Q; was a non- 
linear function of C; alone. On the basis of certain general properties of QO: 
it was shown that for a spherically symmetric system there exists a steady- 
state solution of (1) such that C, = 0 everywhere and tends to zero at 
infinity. This solution represents the spatial accumulation of the solute 
and this phenomenon does not depend upon-the presence of membranes 
or other diffusion barriers. Such a spontaneous accumulation of a solute 
and the resultant “‘structure’’ of the concentration distribution is of some 
interest from the biological point of view. However, certain points in the 
mathematical argument employed by Rashevsky merit comment. 

It is one of the purposes of this paper to furnish such comment. The 
formulations presented here give, under quite broad conditions on Q,, 
some general properties of the steady-state solution of (1). Moreover it 
will be shown that certain conclusions previously drawn (Hearon, Joc. cit.) 
on the basis of constant Q, and constant diffusion coefficients are valid 
under much more general conditions. Finally, certain results of this paper 
will be made the basis for a critique, to be published later, of some ap- 
proximate methods of solving (1). 

1. Simple F ick’s law. In this section it will be assumed that the simple 
Fick’s law, J; = —DiVCi, applies. The diffusion coefficient D, is not 
constant unless so specified. It is further assumed that 7, is a continuous 
vector. From what follows it can be seen that this assumption necessitates 
the additional condition that rQ; is bounded at r = 0. It will be assumed 
that Q; is bounded for 0 < r < ro, where ro is the radius of the cell. These 
conditions and assumptions are referred to as the general conditions. The 
following propositions will now be proved. 


Proposition I. If, in addition to the general conditions, QO, is bounded 
at r = 0, then 


a 
Jp=0s5 r=0, (2) 
and if 
Lim D,= Di¥0 (3) 
then 
Gy 
TAP adr °° (4) 
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Proposition IT. If, in addition to the conditions of J, aD,/dr is bounded 
at r = 0 then 


Lim — =. (S) 


where Q} is the value of Q; at r = 0. 
Proposition IIT. lf I and JI are true then evidently (1/r) dC,/dr is 
bounded at r = 0 and 


Lim 2 2Ce — Or. 


for dr 3.0 oo 


Proposition IV. If, in addition to the general conditions, Q, has a 
constant sign* and D, > 0 for 0 < r < r then C;(r) is monotone and, if 


dC, 


—— =(0 at r=R, (7) 


then Q; changes signs at least once on the interval 0 <7 < sR. 

, Proofs. Under the conditions of spherical symmetry of p= nd (7), where 
n is a unit vector parallel to r, and (1) reads 
1 
re 


V-I,= “T(r =O. (8) 


The solution of (8), under the general conditions, is 


Ji (1) =, f rddr. (9) 

0 

From (9) , 
Lim J.(7) =5 Lim (r Qe) , (10) 


and the truth of J is evident. If (8) is written in full it becomes 


aC; 2 dC; Ox dCi, d (14% 


ie aye Doar a, 


According to (3) and the assumption that Q; is bounded at r = 0, the 
limit, as r — 0, of the first term on the right in (11) exists and is — Qf/Dj. 
From (3), (4), and the boundedness of dD;,/dr at r = 0, the limit of the 
second term on the right is zero. The limit of the left-hand side of (11) is 


* The quantity Q, may necessarily have a fixed sign [e.g., Qi = BCx, 6 constant or n= 
BCx/¢(Cx), @ a polynomial with positive coefficients] in consequence of C; = 0. This quantity 
QO, may potentially change sign, as Qz = a-bCx, a, 6 > 0 and constant (Reid, 1951). In this 
particular case C; is always such that Q; has a given sign. 
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3(d2C;,/dr*). This proves JI and hence IJ. The truth of ZV is evident 
from (9) and the assumption of the simple Fick’s law. 

The conditions required for the proof of J—JV are certain to be satisfied 
in any real physical case. It is to be noted that the conditions D; ~ 0, 
D, > 0in J and JV are thermodynamic requirements. 

In the case discussed by Rashevsky (loc. cit), (4) was imposed as one 
of the boundary conditions and (5) was assumed on the basis of an ap- 
proximate argument involving very particular properties of Q;. It has 
been shown above that under quite general conditions (4) is inherent in 
(8). This redundancy does not alter the conclusions drawn by Rashevsky. 
It has further been shown that, under very reasonable conditions on D; 
and its gradient, (5) is true if Q; is bounded at r = 0. It is also seen from 
IV that if Q;, > 0 everywhere, then dC;/dr < 0 everywhere, as Rashevsky 
assumed. It is clear from (11) that d?C,/dr? does not, contrary to the state- 
ment of Rashevsky, have a fixed sign. 

For future reference, the integral of (9) will be considered here. If it 
is assumed that D, is constant then the integral of (9) is 


Cy (7) =C, (0) <p Gel), (12) 


where 


GOVE [lf hod = Ree 
0 0 0 


Equation (12) follows at once from (9) if it is noted that integration by 
parts gives 


ie [fsa drJan= fr) [sas]ex. (14) 


It is understood that Q;(r) denotes Q;(Ci(r), Co(r),..., 7). It will be 
noted from (12) that if Q, > 0, for all 7, then C,(0) — C,(r) > 0, for 
all r, and conversely if these inequalities are reversed. This result, obvi- 
ously true if JV is true, does not depend upon the assumption, in (12), 
that D, is constant, but is a consequence of the simple Fick’s law. 

2. The general Fick’s law. It is not difficult to establish the counterparts 
of propositions 7-ITTJ when, in the general conditions of Section 1, the 
simple Fick’s law is replaced by the generalization (Onsager, 1945): 


Jum — Do DiVCy © REA, 2y cue 8 (15) 
i 


es 
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If the general conditions of Section 1 are maintained with the exception of 
(15), the following propositions are easily proved. 
Proposition V. If, in addition to the general conditions, every O; is 


bounded* at r = 0 then every VF = 0 at r = 0 and further if each Lim 


. - r>0 
Dx; = D,j exists and the determinant |D,°| + 0 then for all j 
sped epee 
pipe (16) 


Proposition VI. If, in addition to the condition of V, every dD,,;/dr is 
bounded at r = 0 then 


ee Ta 3d (17) 


where 7}; are the elements of the matrix [D,,]~ evaluated at the origin. 

Proofs. Clearly (9) does not depend upon the simple or general form for 
7,(r) and hence the truth of the first statement in V is evident. From 
(10), (15), and the conditions of V 


> di: Lim == 0 PISS, (18) 


from which (16) follows if |D,°| # 0. With (15), (8) may be written as 


D> Disv°Ci+ DY) VDi5-VCi= —Q- (19) 
7 ; 


J 


In view of (16) and the conditions of VJ, (19) gives 


Lim V°C. = — We T3035 (20) 
7 


and, if the form of y?C; be noted from (11), (17) follows at once from (16) 
and (20). The analogue of /// is sufficiently obvious. 
From (9) and (15) 


dr _ SrPti (" p0.d1 (21) 
=a J 70; ’ 


which shows that the constancy of the sign of every Q; is not enough, in 
general, to insure that C;(7) is monotone. In general, then, the counter- 
part of /V is not true. It is to be noted that the condition on |D,;| in V 

* Note here that, in J, if a given Qx is bounded at r = 0 the derivative of the corresponding 


concentration C; vanishes at r = 0. In V, however, if a given Qk is bounded at r = 0 then the 
linear combination (18) vanishes for that k. 
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is included in the thermodynamic requirement that the matrix [Dj] be 
positive definite. 

The following situation has been previously (Hearon 1950a, b) referred 
to: If every Q, is constant then for the simple or general law of diffusion 


— 


Tk (1) = On, (22) 


as is immediately obvious from (9). Thus the direction of the . flux of the 
kth substance is determined by the sign of Q;. In particular 3 , vanishes 
with Q,. For the simple case 


7 Ge 
VO, = D. 3’ (23) 
and for the general case 
r 
VC.= — YT 3: (24) 


Thus the direction of VC; is also determined by the sign of Q:, im the 

simple case. In the general case the direction of VC; is not necessarily that 
> 

of —J;, and is determined not by the sign of Q, alone but by the sign of 


SS T;.;0;. We wish to show here that similar conclusions are valid when 
i 


the Q; are any functions of the C; and r subject to the general conditions of 
Section 1. From (9) 


Ji (1) =f r°O,d r (25) 
0 


and from this, or from (21), it follows that 
r ¥, 
VC = a raf rOjdr. (26) 


Equations (25) and (26) obviously support these conclusions: If Q; has a 
constant sign in 0 <r S ry the direction of the flux J; is determined by 
that sign and, as is intuitively obvious, the flow is outward from the 
cell if Q; > 0 and conversely. The flux vanishes if Q, = 0 for all 7, and 
an inert metabolite exhibits no net flux. If the average rate of metabolism 
of a given solute vanishes, that metabolite does not flow from the cell 
for, from (25), Jk(ro) = 70/3 where Q, is the average metabolic rate: 


3 fr Qxdr/ ce 
0 
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Even if every Q; has a constant sign in 0 < r < rp the direction of VCi 
is not determined by the sign of Q, except in the simple case when 
Tn3 = 0, k #7, Tux = 1/Dy. In fact if the signs of Q, and =T;,0; 
are different the kth species flows against its concentration gradient. Further, 
VC, does not vanish with Q, and an inert solute will exhibit a nonzero 
gradient.* It is easily seen that for the general case (12) becomes, if every 
D,; is constant, 


Cur) =Ce(0) —2 TG s(n) , (27) 


from which it is seen that the sign of the difference C.(0) — C;(ro) is 
determined not by the sign of G;(7o), as in the simple case, but by the 
sign of 27.,G;(ro). This implies that the concentration of a produced metab- 
olite 1s not necessarily higher at the cell center than at the cell membrane 
and that an inert metabolite may be accumulated in or excluded from the 
cell in the sense previously discussed (Hearon, 1950b). 

The above conclusions require no restrictions upon the Q; except those 
noted and, except for (27), no restrictions on the D,; other than those 
required by thermodynamics. This particular demonstration depends 
upon spherical symmetry of the system but it does not seem possible that 
the physical content of the conclusion is so dependent upon the geometry 
of the system. 

This work was aided by a grant from the Dr. Wallace C. and Clara A. 
Abbott Memorial Fund of the University of Chicago. 
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* Stated in another way, if the simple law of diffusion is assumed and Qe = then in the 
steady state C; is a harmonic function and may not exhibit a maximum or minimum value in 
0 <r <7. If (15) is assumed and Q; = 0, then in the steady state C; is not harmonic and, 
according to (16), exhibits an extremum at r = 0 provided every Q; is bounded there. 
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ERRATUM 
To the paper by John Z. Hearon entitled “Comments on the Approximate Solution 
of the Diffusion Equation: I.’ (Vol. 15, pp. 23-31). 
In the above paper, equation (16) on page 27 should read 
1 1 
ir? ee A 


0 t0 


C(O) =C°+ 


fC rodr+24 
0 To 


and equation (21) on page 28 should read 
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The approximation method of N. Rashevsky is discussed and reviewed. It is shown that 
in addition to the explicit assumptions and approximations there is involved the assumption 
that the rate of metabolism is the same at every point in the cell and that the average rate of 
metabolism is different from zero. An expression is given for the error in the approximate 
method when the rate of metabolism is any function of the concentration. It is also shown that 
a solution in the form of that obtained by the approximate method is not possible if the general- 
ized laws of diffusion are assumed to apply. 


1. Introduction. It is well known that the exact solution of the dif- 

fusion equation 
¥4 dC 

—VeJ ar Q = OF ’ (1) 

_ where ie Q, and C are the flux, rate of production (per unit volume), and 

concentration of the solute in question, can be obtained only if assump- 

tions regarding the geometry of the system and the functional form of Q 

are made which often, and for many purposes, render the problem trivial 

from the physical or biological point of view. Attempts to cope with this 

situation have resulted in an approximation method (Rashevsky, 1940, 

1948) which reduces (1) to an ordinary differential equation. The “inte- 


grated”’ form of (1) is 


5 fcav = foav—f Feds, (2) 


where V is the volume of an arbitrary region, S its bounding surface, and 
de a vector element of area with direction of the outward normal to S. 


* A portion of this work was performed while the author was a research participant, Oak 
Ridge Institute of Nuclear Studies, assigned to the Mathematics Panel, Oak Ridge National 


Laboratory. 
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In particular we may choose V and S as the cell volume and the cell 
surface and (2) may be written, provided V is constant in time,* as 


ee 1 > = 
iy ae d 3 
= =) ae do, (3) 


where bars denote averages taken over the volume of the cell, ie., for 
any f, 


set 
faz f sav. (4) 


The final result of this approximation method is equivalent to the asser- 
tion that 


(S) 


where C° is the value which the concentration outside the cell approach- 
es at infinity, A the so-called total diffusion resistance of the cell,f (Ra- 
shevsky, 1940, p. 14, 1948, p. 22), and ~ denotes approximate equal- 
ity. In what follows we will, in accordance with established practice in the 
literature, refer to (3) with (5) as the approximation or the standard approxi- 
mation method. The derivation, discussed briefly below, breaks down in 
certain special cases (see, for example, the discussion of Case iii, in section 
6 of this paper). The relations between the approximation method and 
certain other approximate solutions of (1) are discussed in a forthcoming 
paper (Hearon, 1953b). 

In essence the method consists of approximating the cell geometrically 
by a spheroid of polar and equatorial radii 7; and r2 respectively, comput- 
ing approximate components of I and replacing the surface integral in 
(5) by multiplication of these components by appropriate areas. In com- 
puting the approximate components of Jan assumption is made regarding 
the point at which C = C. This assumption is such that the “average 


*It is to be noted here that the time rate of change of the total amount of solute is given 
by d(VC)/dt = V dC/dt only if V is constant. In the standard approximation method of 
Rashevsky (1940, 1948), V is assumed constant, while in subsequent application of the method 
to cell growth V is allowed to vary with time. 


} This term seems first to have been employed by H. D. Landahl (1939) who neglected the 
contribution of external diffusion (see eqs. [14], [22]). The term has been retained (Rashevsky, 


1940; Rashevsky and Landahl, 1948) even though A, equation (22), is not a property of the 
cell but of the system. 
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concentration” appearing in the J-components* is not the volume aver- 
age, C, appearing in (3). This assumption has been removed by A. S. 
Householder (1942) who treated the “triaxial cell.’”’ The componeiits of 
VC at the membrane in the medium outside the cell are approximated 
by the difference between C® and the concentration at the membrane, 
divided by a quantity 5. This quantity, assumed to be identical for all 
directions, is taken to be the distance at which the external concentration 
would reach C° if VC had everywhere the value it has at the membrane. 
This assumption of the equality of 6 in the three principal directions is 
dropped in the Householder treatment of the triaxial cell. On comparing 
the exact steady-state solution for a spherical cell with constant Q with 
the result of the approximation method, Rashevsky (1940, 1948) assigns 
the quantity 6 the value ro, the radius of the cell. The basis is that, al- 
though this choice is not numerically the “best” choice,t it gives the 
form for A required by the exact solution. For the nonspherical case the 
value assigned to 6 is the smaller of 7; and r. (Rashevsky and Landahl, 
1940). 

2. The problem. Although it can be shown for the spherical case with 
constant Q that the final result is relatively insensitive to the various geo- 
metric approximations and to the choice of the value of 7 at which it is 
assumed that C(r) = C, in general the error propagated into the final 
result is not obvious. What does appear obvious is this: In general, A 
will depend upon the functional form of Q. This is easily seen by compar- 
ing the exact steady-state solution for C for some simple cases with the 
result 

C=C°+A0 (6) 


obtained from (5) and (3), or by considering the explicit formulation 
given by A. Weinberg (1939) for the case 


O. = AyrCi + A2C2 , k=1,2. (7) 


Moreover, it is physically obvious that the actual value of the quantity 6, 
discussed above, will depend upon Q. 
3. Purpose and scope of paper. Since comparisons of the approximate 


* The ‘‘average concentration” used by Rashevsky (1940, 1948) is actually [C(0) — C(ro)]/2 
for the spherical case. This value of C(r) occurs at r = rol /2, if Q is constant. The ve of C 
occurs at r = 79x/6/+/10 and the radial average at ro/+/3. The ‘‘mean vC-components com- 
puted from these ‘‘mean concentrations”’ as [C* —-C(ro)]/(ro/p), where C* is [C (0) — C(ro)]/2, 
C or the radial average and ro/p is the appropriate distance in each case, do not differ greatly. 
It does not appear likely that such agreement will be secured when Q is not constant. 


+ Evidently exact agreement is obtained if 6 = D,/2h + (3/2 — D./5Di)ro. 
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and exact solutions have been made (Rashevsky, 1940, 1948, the spherical 
case; Young, 1939, the ellipsoidal case) only for constant O, the effect of 
the functional form of Q is worth considering. 
It is the purpose of this paper to show that for the steady-state spher- 
ical cell 
C=C4+(A+e00, O#0 (8) 
=C—770/6D;, Q=0, 
where A is the total diffusion resistance for the spherical cell with constant 
QO (eq. [22] below), D; the diffusion coefficient inside the cell, and ¢€ a 
“correction term.’’ The sign and magnitude of e depend upon the func- 


tional form of Q: 


Ps 
ro 


1 OP 


€ (1— PO/ Q) ; 
in particular, of course, e = 0 if Q is constant. 

It is not the purpose of this paper but of a later one (Hearon, 1953b) 
to compute ¢ for any given case. It may be possible to give some general 
conditions under which the standard approximation method is valid. In 
any event the present formulation exhibits in a single additive term the 
discrepancy between the exact and approximate solutions. The assump- 
tion that Q is of such a form (or involves parameters of such a magnitude) 
that ¢ is negligible with respect to A must be added to the other assump- 
tions in the approximation method. Previously no explicit restrictions 
have been placed upon Q. Fortunately, however, in many applications 
of the approximation method it has been assumed that Q is constant. 

4. The “integral equation” solution. It has been shown elsewhere 
(Hearon, 1953a) that the solution of (1) for the steady-state spherical 
case is 


Jet rOdr. (9) 
This solution is quite independent of the functional form of Js viz., does 
not depend upon the simple or generalized (Onsager, 1945) form of 
Fick’s law and is independent of the presence of convection and force 
fields. If the simple Fick’s law and the absence of convective flow and 
force fields are assumed then, provided that the diffusion coefficient, D;, 


is constant, the solution of (9) for the concentration, C(r), is (Hearon 
1953a) 


C(r) =C(0) Se), (10) 
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where 
G AY Pog ee fe 
(r) i (r—2) O(A) dd eal AQ (A) ddd)y. (11) 


In (9), (10), and (11) Q may be any function of C andr provided (Hearon, 
1953a) it is continuous for 0 <r< ro and rQ is bounded at r = 0. In 
fact Q may depend upon concentrations of other substance and, for 
simplicity, we denote by Q(r) the function Q(C1, Co, ... ny Todt Gale) 
is the concentration in the medium, then, since it is assumed that Q0=0 
everywhere outside the cell, 


C.(r) =C°+ b/r, b constant . (12) 
The boundary conditions are 
C(r) bounded at r = a 


Pele o Oh Mel pes (13) 


and 
> o> hes 
neJ=—D, ie 
at r= 1p i) 
=h(C—C,) 


where 7 is a unit vector having the direction of r, D, the diffusion co- 
efficient in the medium, and / the permeability coefficient of the mem- 
brane. 

Conditions (13) are included in (10) and (12) (cf. Hearon, 1953a), and 
C(0) and 6 are arbitrary constants to be fixed by conditions (14). From (9) 
and (12) the first of conditions (14) gives b/7ro as 


il To 
a ee 5 
b/ ro Da, tp rOdr. (CUS) 
With (10) and (12) the second of conditions (14) evaluates C(0) as 


1 7 b il 
C (0) SCs ae) rQdr ++ 5G (1) . (16) 


Equations (15) and (16) give the arbitrary constants in (10) and (12) as 
determined by the conditions (14). Unless Q is a function of r alone, (10), 
with C(0) from (16), is an integral equation for C(r) of an extremely 
complex type. ae | 

5. Comparison to approximation method. For our purposes It 1s con- 
venient to express the constants b/ro, and C(0), as well as certain quanti- 
ties which appear in the expression for C, in terms of Q, of the volume 
averages of certain functions of Q, and of the diffusion properties of the 
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system. This is a lengthy but simple procedure which will now be out- 


lined. 
If the system is spherically symmetric (4) becomes 


faa f cae 


It is then found, after noting (15) and the frst form of G(r) in (11), that 


2 3 2 
ao 2G Oe) ee 
C (0) =O +37 Ot an Ut 3p eee ee iad 
From (10), it follows that 
as 3 To 
C=C) aah Glridr. (18) 


Denote by F (ro) the integral in (18); then after two integrations by parts, 
in a manner which will be clear from (19), there results 


d To 2 
Fre) =26(n) — ee waa ent dr. (19) 


If G(ro) in the first term of (19) is replaced by the first form from (11) 
and G(r) in the second and third terms by the second form from (11), it is 
found that 


Age 


Det, se 


3 fee es vere. <8 
~ Dir rH Oe, Ore ap FO). (20) 


Equations (17), (18), and (20) give 
2 2 
ex 10. sahtare et en) tical aces 
Gmc: (eee) Oey nn: a 


There are now three cases to be considered: ne 
(i) If Q is constant, (21) and the relation 72 = 3r6/5 immediately 
yield (6) 
=C? + AO ’ 
where 


2 
9 


15D;° 


2 
‘9 ro 


h (22) 


Equations (6) and (22) constitute the exact solution for @ with which, 
except for numerical factors in (22), the approximation method agrees. 
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(ii) Lf Q is not constant and OQ + 0, (21) can be written as 


C=C+ (A+600, (23) 


which is the first of equations (8) in which A is given by (22) and the 
correction term e is 


2 
9 


aap COLT 0! . (24) 


; = If Q = 0, then the second of equations (8) results directly from 
21). 

6. Discussion. While it is not the purpose here to explore in detail the 
limits which can be set upon the correction term e, a few obvious points 
will be noted. Evidently if Q is a prescribed function of r alone ¢ can be 
computed at once. This will seldom be the case in practice (cf. discussion 
in Hearon, 1953a) but some idea of the values of e, when Q departs in a 
given manner from constancy, can be gained. Since r is inherently positive 
a standard theorem gives 


(1—Qu/Q) < (1— 770/72 Q) < (1-O,/0), (25) 


whence Qy and Q,, are the maximum and minimum values taken on by 
Qin JV. 

Case (iii) above is admittedly a very special situation; however, the 
standard approximation method yields, in this case, a completely er- 
roneous result except for the trivial case Q = Q = 0. It is necessary to 
begin anew and derive an approximate solution in this case as well as 
numerous others (see treatment of cell polarity, Rashevsky, 1948). In 
order that Q = 0 the substance must be produced in some regions of the 
cell and consumed in others. A formally identical situation has been dis- 
cussed by H. D. Landahl, 1943 in which it was assumed that k = 0. From 
(9) and (14) this assumption clearly requires 0 = 0. Moreover, although 
we shall not discuss it further, this would seem to be an important case. 
For it is well known that many important metabolites do not flow ap- 
preciably, if at all, from the cell [note from (9) that if Q = 0 then J = 0 
at r = ro; cf. Hearon, 1953a]. An important defect of the original dif- 
fusion force theory of cell division (Rashevsky, 1948) is that a net outflow 
of material from the cell is required for division to occur.* It was this 
consideration which prompted the investigation of Landahl (1943). 

* It can be shown (Hearon, 1950, 1953a) that if the generalized law of diffusion is assumed 
it is not necessarily true that C(0) > C(ro) if Q > 0 and conversely. It is clear from this paper 
(eq. [27], etc.) and earlier work (Hearon, 1950, 1952) that the general laws of diffusion, which 


are certainly qualitatively important where electrolytes are involved, require considerable 
revision of the diffusion-force ideas of cell division as originally presented. 
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The demands which different applications make upon the accuracy of 
the approximation method will vary considerably. In this connection it 
is to be noted that in practice (Rashevsky, 1948; Hearon, 1949) when it 
is desired to formulate O as a function of C° and the diffusion properties 
of the system, Q is expressed as a function of C by assuming that 


Q(C) =O). (26) 


Clearly (26) involves a rather gross approximation unless Q(C) is a linear 
function of C. If (26) is to be accepted one cannot logically require very 
high accuracy, in other respects, from approximation methods. It is clear 
that whenever the magnitudes of diffusion and osmotic forces are com- 
puted* by the approximation method the results will be in error except 
when (Q is constant. 

It should not be difficult to extend the present formulation to other 
systems of simple geometry, for example, the cylinder with assumed axial 
symmetry. In particular, the “metabolizing slab,’”’ approximated in prac- 
tice in the usual tissue slice experiment, would be of interest. A compari- 
son of the so-called limiting tissue thickness for tissue slices as computed 
by the exact and approximate methods has been made (Rashevsky and 
Landahl, 1940) for constant Q. Differences between the predicted and ob- 
served thickness are seriously discussed (Shorr, 1943) despite the fact that 
the assumption of constant Q is manifestly too crude an approximation 
for small C°. 

It is simple to show (cf. eq. [27], Hearon, 1953a) that extension of the 
present formulation to the case of Onsager’s generalization of Fick’s law 
(Onsager, 1945) gives, instead of (21), 


Ce=C2+ (Ret Ox) Or » (27) 
where 


To pe r 87) ro tH 1 

R,= D oF AiQit> TiQit Ts Ti0;) Or 
(28) 

Tool Sex a iar 

= 19 Oye Pe — 7°Q;/ 7 Qj) 


[Ti] = (Dil, (Tj) = (DRI, — [Hl = [hal 1 (29) 


* There is no advantage to an actual computation in the present case for the sphere is 
necessarily a stable configuration, but it is plain that if C and C(ro) are in error the diffusion 
and osmotic forces are correspondingly in error. If C(r) is not monotone the approximate 
formulation may even yield the wrong sign (direction) for these forces, 
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and the Di;, Di;, and Ay; are the internal and external diffusion coef- 
ficients and the “generalized permeability coefficients” (Hearon, 1950). 
When every Q; is constant 6; vanishes. If [Dj,] is the diagonal matrix 
[5:; Di], 9 reduces to ¢, as given by (24). Clearly, however, im general 
there exists no counterpart of the approximation method in the sense that no 
equation of the form (6) or (8) exists where A [replaced by R; in (27)] 
depends only upon the dimensions and diffusion properties of the cell. In 
the special case of only two substances, it is necessarily true that 0; = 
—Qs2 and it is evident that under these conditions Ri, k = 1, 2, depends 
only upon the diffusion properties of the system. 

The above statement does not mean, of course, that no approximate 
solution based upon the kinds of assumptions found in the standard ap- 
proximation method exists. The difference in form between (27) and (8) 
is worth noting however. Moreover, it would be of interest to study (27) 
further in order to determine the importance of the “‘correction terms,” 
6,, when coupling by diffusion (Hearon, 1950b) can occur. 

This work was aided in part by a grant from the Dr. Wallace C. and 
Clara A. Abbott Memorial Fund of the University of Chicago. 
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Rate equations for the enzymatic oxidation of succinic acid are derived on the assumption 
that when a single molecule of substrate combines with an enzyme molecule, it can do so with 
either one or two sites on the enzyme, and that oxidation occurs only in the second case. In 
addition it is assumed that the product of the reaction, fumaric acid, combines reversibly with 
the enzyme. With certain enzyme preparations the data fitted such an equation satisfactorily. 
In others the rate was that of a first-order reaction, but addition of cytochrome changed it to 
the former type. It was concluded that the transfer of hydrogen to oxygen was a first-order 
reaction and dominated the whole rate when enzyme preparations were used which had been 
washed relatively free of cytochrome. When the limiting factor was succino-dehydrogenase the 
rates followed the new equation. Criteria for recognizing noncompetitive inhibition are given, 
and inhibition by di-tertiary butyl peroxide was shown to be of this type. 


The oxidation of succinic acid in the presence of tissue succino- 
dehydrogenase and cytochrome was found to be a reaction of the first 
order by E. Stotz and A. B. Hastings (1937) who used an enzyme prepara- 
tion which had been carefully washed so as to free it as far as possible from 
fumarase. A similar conclusion had previously been reached by F. P. 
Mazza and T. Laurenza (1934), but their data frequently show a definite 
trend in the values calculated for the velocity “constant,”’ and are there- 
fore less convincing. Before this, however, other workers (see, e.g., Ahl- 
gren, 1928) had measured the rate of the reaction by the time of decolor- 
ization of a standard amount of methylene blue in vacuo. This is justifiable 
only if the rate of the reaction is approximately constant throughout, 
which is not true for reactions of the first order where the rate decreases 
logarithmically. No direct proof of the constancy of the rate was available 
for tissue enzymes, although it has been shown for the oxidation of 
succinic acid by resting bacteria (Yudkin, 1934; Quastel, 1926) and for the 
oxidation of xanthine by an enzyme system obtained from milk (Dixon 
and Thurlow, 1924). Moreover most workers, beginning with F. Batteli 
and L. Stern (1911), agree that high concentrations of substrate inhibit, 
an effect difficult to explain on the basis of a first-order reaction. 


33 


34 G. S. EADIE AND F. BERNHEIM 


Our own data in certain experiments indicated rates differing markedly 
from those of a first-order reaction, and it was therefore necessary to 
consider another hypothesis. The symmetry of the succinic acid molecule, 
and the inhibition by high substrate concentrations suggested that dehy- 
drogenation might occur only if the succinic acid molecule became at- 
tached to the enzyme at two points. Such a suggestion had been made by 
Mazza and Laurenza (loc. cit.) although they did not pursue its kinetic 
consequences. In addition, the usual assumptions are made that the com- 
bination of substrate and enzyme is a reversible reaction reaching equilib- 
rium very quickly as compared with the breakdown of the compound to 
yield the products of the reaction. Because the reaction is inhibited by its 
product, fumaric acid, it was also assumed that fumaric acid could com- 
bine with the enzyme in exactly the same ways as succinic, but with dif- 
ferent dissociation constants. These reactions would give rise to a series 
of simultaneous equilibria between enzyme, substrate, and product. The 
possible equilibria involved may be represented by the following equa- 


tions where HSuc_ stands for succinic acid, HFum_ for fumaric acid, 


and E for the enzyme system. The bonds represent possible points of 
attachment: 


wee PA eka —-E (1) 


The product of this reaction results from attachment at one site; but at- 
tachment at the second site may also occur, giving rise to another equilib- 


| ne x 


HSuc-E = HSuc=E, (2) 
or the product of reaction (1) may combine with a second molecule of suc- 
cinic acid 


NOR ca i, 


HSuc—E + #£HSuc HSuc—E—HSuc , (3) 


or with a molecule of fumaric acid, 


\ Se ae if 


HSuc—E + ‘HFum = HSuc—E + HFam* : (4) 
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Fumaric acid may be substituted everywhere for succinic acid, giving us 
three more equations: 


Sad a A = \ a 


= HFum-E , (oS 


HFum—E = HFum=E, (6) 


< Sa 


HFum—-E + HFum 2 None ee (7) 


It was assumed that these equilibria become established almost in- 
stantaneously after mixing and that the rate of dehydrogenation is pro- 
portional to the concentration of the product of equation (2). The reverse 
reaction, the reduction of fumaric acid to succinic, might, similarly, be 
proportional to the concentration of the product of equation (6), but 
under the conditions of these experiments, v7z., 20% oxygen, reduction 
does not occur, and reverse reactions may therefore be neglected. 

To express these equilibria in quantitative terms let the concentration 
of the total unoxidized succinic acid be [A]; the concentration of free suc- 
cinic acid is theoretically less than this, but under most experimental 
conditions, including those of the experiments reported below, the differ- 
ence is very small, and may be neglected, as in the Michaelis-Menten 
theory. Let the free enzyme concentration be [Z,] and the concentration 
of the product of reaction (1) be [C]. Then 


[A] + (E)] =2IC1, (8) 


where 1/e is the dissociation constant of the compound C. Similarly, 
let [D] be the concentration of the product of reaction (2) and 1/y its dis- 
sociation constant; then 


Te eae Bale (9) 
“ig 


Let [F] be the concentration of the product of reaction (3) and 1/6 its dis- 
sociation constant; then 


[C]-[A] =<IF]. (10) 


Let [A’] be the concentration of total fumaric acid, assumed to be prac- 
tically equal to that of the free acid, [G] the concentration of the product 
of reaction (4), and 1/¢ its dissociation constant; then 


(C] «(4/1 =2161 . (11) 


36 G. S. EADIE AND F. BERNHEIM 


Similarly, we may write for reactions (5), (6), and (7) 


[A’] + [Ey] =5IC'l, (12) 
il 

(oc) = (D1, (13) 
7’ 

[C’] + [A oe (14) 


From (10) and (14), putting 6’/6 = 


a ic [A’] z 
Similarly, if y’/y = wand e’/e = , 
(oy 
LP ear (16) 
and 
ee [A’] i 
len= lela dr. (17) 
Substituting for [C’]/[C] from (17) into (16) and (15) we find 
[D'| = (D144, (18) 
ee LAA 
LP’). = Fl aq ATP: (19) 


The free enzyme concentration at equilibrium must equal the total en- 
zyme concentration less the concentrations of the various compounds 
containing it, and hence, because of equation (8), 


[Eyl =< (B). = CIC] 7 (Geet Dae 
CUR hens WT) =. 1G) ee oe eal ee 


Substituting for [C’], [D’], [F’], and [G] from (17), (18), (19), and (11), 
then eliminating [C] and [F] by the use of (15) and (16) we — an ex- 
pression, which, solved for [D], becomes 


2 ie 
[A]+ (1+ 7+6[A])+(0+yrAu+ddr7[A’]) + 


(20) 


p= 


mCUE I ee 


spayed 


But the velocity is assumed proportional to [D]. Therefore if « = [A’], 


dx 
a Fl?) . (22) 
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Letting a = [A] + [A’] and integrating (22) with its initial conditions 
x = 0, t = 0 we obtain 


a 


a, Px TMGgsesy +Rx? —-1, (23) 
[E] ykP = 1+ y+ 6a—2X(1+yu4 ba), (24) 
LE] yRQ=e 1+ ad(1+ Au ba), (Zo) 
[E] yRR=3(£—6—6»). (26) 


The initial velocity in solutions containing no preformed fumaric acid 
can be obtained by equating [A’] to zero in (21), multiplying [D] thus 
obtained by &, and rearranging: 


ky(A] [E] 


"TE Ut 7+ 64] (AT 2b 


By differentiating this with respect to [A] and equating to zero it can be 
shown that a maximum velocity occurs when 
ie 
[A]o= ae: (28) 
€ 

Equations identical in form but with slightly different constants can be 
obtained if it is assumed, for example, that fumaric acid does not combine 
at all with the enzyme, or that it combines at one site only, i.e., that the 
equilibrium represented in equation (6) does not occur. Another alternate 
set of hypotheses is that oxidation of succinic acid occurs if the substrate 
molecule is attached to the enzyme at one site only, provided also that 
the neighboring site is unoccupied. This also leads to an equation identical 
in form, but again with slightly different constants. 

Noncompetitive inhibition will result in new values for the constants 
of equation (23), P’, Q’, and R’, such that P’: P = Q':0 = R’:R. This 
follows from equations (24) to (26). In general the net result of noncom- 
petitive inhibition is to reduce the concentration of active enzyme (or the 
value of Vmax Which is proportional to it) to a value equal to the original 
value divided by (1 +72/K,). Since, for example, [Z]ykP is constant 
[eq. (24)] and since y and & cannot alter, when [E] is divided by this fac- 
tor, P must be multiplied by the same factor to keep the value of the 
whole expression constant. Similarly, Q and R must be multiplied by the 
same factor. This is not true of competitive inhibition in which the disso- 
ciation constant of the enzyme-substrate compound is apparently altered. 
This would affect the constants P, Q, and R unequally. 
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If the values of the constants are such that 


¢=6(14+»), (29) 
then R [eq. (26)] becomes zero, and equation (23) reduces to 
Pet Oly (30) 
a—xX 


which is the same as the equation of L. Michaelis and M. L. Menten 
(1913) or D. D. Van Slyke (1942). If, in addition, P also becomes equal to 
zero, the equation reduces to that for a first-order reaction, but this im- 


plies that 
1+y+6a=dA(1+ yu) + Adve, (31) 


and the presence of a, the initial substrate concentration, in two of the 
terms makes it unlikely that this condition could often be fulfilled. 

Rat liver was used as the source of the succinoxidase. Five gm. of 
homogenate was suspended in 50 ml. of 0.025 M K-Na phosphate buffer 
at pH 7.8 and centrifuged. The solid was resuspended and centrifuged 
again, and this was repeated twice. It was then taken up in 0.05 M buffer 
pH 7.8 and placed in Warburg flasks with a 2.0 ml. fluid volume. The oxy- 
gen uptake in the presence of succinate added from the side arms was 
measured in the usual way. 


RESULTS 


When the data from a number of experiments were plotted, it was seen 
that only a small number of curves were logarithmic in shape as of a first- 
order reaction (Fig. 1). Others differed in rising more steeply and remain- 
ing linear until some three-fourths or more of the substrate was oxidized; 
the rate then fell off very sharply (Fig. 2). The remainder were interme- 
diate in shape between these two. Curves of the second, steeply rising, 
type could be fitted by equation (23) by a suitable choice of constants; the 
smooth line in Figure 2 was drawn from such an equation with constants 
P = 0.061, Q = 12.78, and R = —0.000053. Values of the same order of 
magnitude were found to fit other curves. Unfortunately from these three 
values it is impossible to compute values for the six or seven original dis- 
sociation constants. It was also seen that the shape of the curve depended 
on the enzyme preparation used; this suggested that washing might have 
removed a substance responsible for the difference. From the work of 
Stotz and Hastings this might be fumarase; in the absence of this enzyme — 
fumaric acid would accumulate and slow the reaction by competing with 
the substrate for the enzyme, but not, under the conditions of these ex- 
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Ficure 1. Rate of oxidation of succinic acid. The logarithm of the still unoxidized amount 
of acid, measured as cu. mm. of the equivalent oxygen, is plotted against the time in minutes. 
The straight line is drawn from the equation Log(a — x) = 2.8976 — 0.00854. The reaction 
is therefore of the first order. 
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Ficure 2. Another type of oxidation curve for succinic acid. Oxygen consumption in cu. 
mm. is plotted against time in minutes. The growth curve is drawn from the equaticn 
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Ficure 3. Rate of oxidation of succinic acid. Log(a¢ — x) is plotted against ¢ as in Figure 1. 
The same enzyme preparation as in Figure 1 was used, except that 0.1 mg. cytochrome C was 
added. The points do not fall on a straight line and the reaction is therefore not of the first 


order. 
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Ficure 4. The same data as in Figure 3. Oxygen consumption in cu. mm. is plotted against 
time in minutes. The smooth curve is drawn from the equation 


0.04844 19.88 log —“— — 0.00002 12" = 1. 
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periments, by driving the reaction backward. Experimentally it was 
found that the more completely washed preparations gave first-order 
reaction curves, and it therefore seemed unlikely that the change in 
fumarase content of the preparation could be the explanation. The cii- 
zyme consists not only of succino-dehydrogenase but also of cytochrome 
C, and the latter is more easily washed away than the former. Therefore 
0.1 mg. of cytochrome C was added to a well-washed enzyme preparation. 
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Ficure 5. Effect of di-tertiary butyl peroxide on the velocity. The factor, which is the re- 
ciprocal of [1+(z/Kz)] (see text), is plotted against mg. of inhibitor added to 2 cc. of the re- 
action mixture. 


Data for an experiment with this preparation without added cytochrome 
are given in Figure 2 where log (a — «) is plotted against time; points 
fall on a straight line, and the reaction is therefore of the first order. In 
Figure 3 data are plotted for the same enzyme after addition of cyto- 
chrome; the points no longer fall on a straight line, and the type of reac- 
tion has changed. The data for this experiment are fitted by equation (23) 
as shown in Figure 4. Addition of cytochrome is thus sufficient to change 
the curve from that of a first-order reaction to that demanded by the new 
theory involving double contact of enzyme and substrate. The significance 
of this appears to be as follows. Oxidation of succinic acid proceeds in two 
steps: dehydrogenation catalyzed by succino-dehydrogenase, and com- 
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bination of the hydrogen with oxygen catalyzed by cytochrome. Under 
the conditions of these experiments neither part of the reaction can occur 
(except to a negligible extent) without the other. When the concentration 
of cytochrome has been considerably reduced with respect to that of 
succino-dehydrogenase by washing the enzyme preparation, the oxidation 
of the hydrogen will be relatively slower than dehydrogenation and its 
rate will therefore dominate the rate of the complete reaction. Addition 
of an excess of cytochrome reverses the relative concentrations and de- 
hydrogenation becomes relatively slower, and therefore dominates the 
total rate. This implies further that oxidation of the hydrogen is a first- 
order reaction, and that dehydrogenation of succinic acid involves a 
double contact of enzyme and substrate according to one of the sets of 
hypotheses discussed above. 

Certain fatty acid peroxides inhibit the succinoxidase (Bernheim, 
Wilbur, and Kenaston, 1952). Di-tertiary butyl peroxide was used to deter- 
mine the type of inhibition by this class of compound. It was found to be 
noncompetitive; in every case P, Q, and R of equation (23) were found 
to be altered to the same extent, and the factor by which they had to be 
multiplied varied with the amount of inhibitor added. Figure 5 shows this 
relation; the slope of the line is equal to 1/K;, and in this experiment K; 
was found to be 4.6 X 10-%. 
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This paper deals with the following question: Which distributions of radiosensitivity in a 
population can lead to an exponential survival curve? The problem is solved exactly, with sta- 
tistical fluctuations in dose fully accounted for. It is shown that only an exponential distribu- 
tion of sensitivities can give rise to an exponential survival curve. 


In many types of radiation effects, the number of unaffected organisms 
in an irradiated population decreases exponentially with the dose. It has 
been shown (Wijsman, 1951) that if all the members of the population 
have the same radiosensitivity, such a dose-effect relation implies that the 
effect is always achieved by a single radiation particle (or quantum). But 
one may assume, alternatively, that the effect may be the cumulative ac- 
tion of many hits, and that the exponential survival curveT is due to a 
suitable spread in the sensitivity of the organisms. The purpose of this 
note is two-fold: first, to determine the sensitivity distributions which 
could give an exponential curve, and secondly to consider the possibility 
of deciding which hypothesis applies in any case. 

If the effect is cumulative and requires a great many hits, the organisms 
all receive essentially the same number of hits. A small increment of dose 
then reveals those organisms which are sensitive to the corresponding 
range of hits. Therefore the sensitivity distribution which yields an ex- 
ponential survival curve is given by the derivative of that curve, which is 
also exponential. 

However, the same argument no longer applies when the effect requires 
few hits. In this case, for a given significant dose, the individual organisms 
have received a widely varying number of hits, in accord with statistical 
fluctuations. Therefore an incremental dose no longer critically affects 

* Post-doctoral Research Fellow of the National Cancer Institute. Present address: Physics 
Dept., Memorial Center, New York 21. 

+ Hereafter the terminology is that used for lethal effects, but the discussion applies equally 
well to all radiation effects. 
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and reveals organisms of a narrow band of sensitivities, and the derivative 
loses its significance. The following analysis solves the problem of the 
sensitivity distribution required to yield exponential curves in the gen- 
eral case. 

In this discussion a “hit” on any organism is defined simply as the pas- 
sage of a radiation particle through some geometrical region associated 
with that organism and believed to be basically connected with the effect. 
The hit-region may be a small part of a cell or it may even be a region sur- 
rounding the organism. A “one-hit case”’ refers to an effect due always en- 
tirely to the action of a single particle, although this particle may be 
preceded by totally ineffective hits. The term “cumulative-hit case” in- 
cludes all other situations. 

The following notation is used: 


D = dose, expressed as the total number of particles 
incident on the population, 
p = probability of a hit per particle per organism, 
F(h) = probability that an arbitrarily chosen organism 
survives h hits,* 
P(p,D) = probability that an arbitrarily chosen organism 
survives a dose D. The function P thus repre- 
sents the observed survival curve. 


The probability that an organism receives exactly / hits and survives 
is given by 


es <- ase: 
Pil =f) PS 5 RG with eal 14) EP i 
Therefore 
D 
P(p, D) = >) ph(1— p) 0G, F (h) . (1) 
A=0 


For a one-hit case, the survival function has the form 
PD) Ee oe (2) 
where p’ is the probability of an effective hit per particle per organism.+ 


* The probability F’(h) describes completely the distribution of sensitivity in the popula- 
tion. The fraction of the organisms which succumb on their Ath hit is S(h) = F(h—-1) — F(A). 


t The conventional theoretical curve for the one-hit case, P(D) = e~?’D, is an approxi- 
mation if p’ is given its usual physical significance. The exact general form is (1 — p’)?. 
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In order to determine which functions F(h) may lead to a “one-hit” sur- 
vival curve, we equate the forms (1) and (2): 


D 
P(p, D) = s p*(1 — p) PCF (h) = (1— 9’). (3) 
h=0 


By applying relation (3) to all doses from zero up to some maximum 
dose D, one obtains (D + 1) equations in the (D + 1) unknowns F(A), 
h ranging from zero through D. Thus F(z) is completely and uniquely 
determined. The solution for D = 0 gives F(0); the solution for D = 1 
then yields F(1) and so forth. When carried out this procedure soon sug- 
gests that the general solution of (3) is 


Bh jr ke" with p=1-F. (4) 
Substitution of this function into (1) gives 
D 
P(p, D) = >) (1—p)?-*(p—p') "nC = (1p) + (p—9') 19 = (1-9), 
h=0 


so that this form for F(/) is indeed the solution of (3). 

Equation (4) has this significance: If a population is characterized by a 
distribution of sensitivities F(z) = k’, and by a hit probability p, the 
observed survival curve will be identical with that to be expected for 
“one-hit”? organisms characterized by a probability for effective hits 
pap lh — ky. 

The function F(z) is exponential; F(z) = e*™*. Furthermore, the func- 
tion f(z) is exponential: 


f(h) =F (h—1) —F (h) =h1_b=(G-1)¢ oo 


Thus, if either F or f is used to describe the sensitivity distribution, the 
one and only sensitivity distribution which can yield an exponential survival 
curve is an exponential distribution. 

The form (4) also includes, as it must, the special case of one-hit theory 
with every hit effective. For this case F(/) is given by 


1,h=0 
FO). et. 
This is of the form (4) with & = 0, and it may be noted additionally that 
the substitution of these values for F(#) into equation (1) leads directly 


to the ‘one-hit’? form of equation (2). 
It is instructive to consider the different biological situations that may 
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give an exponential survival curve, and to consider the behavior of the 
function F(/) in each case. 

The simplest category is the one-hit case with every hit effective. This 
occurs, for example, if there is an indispensable molecule which is func- 
tionally destroyed whenever it is hit. In this case, as just noted, F(Z) = k” 
with k equal to zero. 

A second category is the one-hit case where each hit has some proba- 
bility, less than one, of being effective. This is true of many biological 
situations. For example, if the disintegration of a molecule depends on 
ionization density, the fluctuations along the ionization tracks may make 
some hits ineffective. Again, if the action is indirect—activated radicals 
diffusing to and disturbing some critical body—some hits will be ineffec- 
tive because no radicals happen to reach that body. Another possibility 
in this category is that the “‘hit-region”’ is too broadly defined; for ex- 
ample, if every collision with some vital molecule is called a “hit,” while 
only part of the molecule is really sensitive, then the “‘hits”’ outside of this 
part are ineffective. For every case in this second category, we may define 
a probability (1 — &) that any particular hit is effective; the probability 
that the organism survives / hits is then F(h) = k*. 

A final category is the cumulative-hit case with the peculiar sensitivity 
distribution of the form F(z) = k*. This situation is practically indis- 
tinguishable from the one-hit case precisely because F() is defined in this 
case to be the same as the function F(/) which arises naturally, as shown 
directly above, out of the one-hit situations. 

We come now to the question of choosing between the one-hit and the 
cumulative-hit hypotheses when an exponential survival curve has been 
obtained. 

A powerful argument against the cumulative-hit hypothesis is the a pri- 
ori unlikelihood of encountering an exponential spread in sensitivity. This 
artificial distribution has not been found in nature. 

Further important evidence may come from rate-of-dose studies. 
Cumulative effects may depend on rate of dose, while one-hit effects 
cannot. 

Additional evidence may come from experiments with very precisely 
aimed radiation. An exponential survival curve determines just one quan- 
tity: p’ = p(1 — k), which in all cases is the probability per particle per 
organism of consummating the radiation effect. In most experimental ar- 
rangements, p is very small and not well known, so that the obtained 
value of p’ tells nothing about the magnitude of &. But if one was dealing 
with a one-hit case where almost all hits were effective, and could aim the 
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radiation well enough to make # almost unity, the data would show #’ to 
be almost unity and hence & to be much less than one. This would greatly 
reduce the appeal of the cumulative-hit hypothesis, since cumulative ef- 
fects are most often thought of as involving large numbers of hits, which 
correspond to larger values of k. Furthermore, if & is shown to be quite 
small, this means that most of the organisms, at least, succumb to a 
single hit. 

One might hope also to obtain useful evidence by studying fluctuations 
from the average survival curves. But this procedure cannot be helpful. 
The probability of having exactly » survivors in a population of WV or- 
ganisms after a dose D is 


iG SE AU he aod ED of Bala OP 


Therefore, if cumulative and one-hit hypotheses predict the same average 
survival P, they also predict the same fluctuations from that average. 

In conclusion, an exponential survival curve can be due to a cumulative 
effect only if the organisms are characterized by an exponential distribu- 
tion in sensitivity, which seems quite unlikely. Studies with varied dose 
rates and with very precisely aimed beams of radiation can give some 
further basis for choosing between one-hit and cumulative effect hy- 
potheses. 
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An approximation method is introduced which enables a number of diffusion-type problems 
to be solved in an approximate but simple manner. Many cases require only the solution of a 
simple first-order differential equation. The method is applied to a number of cases in which 
the exact solutions are available. A comparison shows that the method is quite satisfactory in 
these cases. The method is applied to diffusion problems with rate of consumption propor- 
tional to concentration or to the square of the concentration. In the latter case, the result ob- 
tained is essentially the same as that found by H. G. Landau (1950) after elaborate calcula- 
tions. 


The approximation method introduced by N. Rashevsky (1948, Ap- 
pendix) for the solution of diffusion problems has proved to be useful for 
the solution of a variety of problems. It is the purpose of the present note 
to extend and improve this method. We shall illustrate the method by 
solving a few simple cases and comparing some of the results with the 
corresponding solutions of the partial differential equations. 


Case 1. One-dimensional, semi-infinite region with no sources c = Co, 
x < 0. Initially c = 0,t = 0,x >0. 

Consider a semi-infinite region containing no sources in which the con- 
centration is maintained at a constant value cy at x = 0, starting at time 
t = 0. Initially the concentration is everywhere zero. 

We make the approximation that at any time #, c(x,t) is given by a 
linear function (Fig. 1) so that the concentration becomes zero at some 
distance r from the origin. This distance 7 is some function of time but is 
zero only for t = 0. Hence we have, for ¢ > 0, 


% 
(x) =a(1-) for TES 


G=0 for se 1A ) 


Expression (1) satisfies the boundary conditions that c(x,0) = 0 for a > 0 
and c(0,f) = co unless r(f) = 0. That the first boundary condition is 
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satisfied can be seen as follows. Since r is zero for t = 0, then for ¢ = 0 the 
statement x > 0 implies that x > 7. Hence, from the second expression 
(1), c is equal to zero. The flow F across the boundary x = 0 is given by 
the product of the diffusion coefficient and the negative of the concen- 
tration gradient at x = 0. Then the time integral of F is the total amount 
which has passed the boundary « = 0. This amount must be equal to 
the total amount of material on the right of x = 0 at the time ¢. In this 


) Xi A(t) 
x——> 


FIGURE 1 


case the total is the product of the average value of ¢ and the distance r, 
or cor/2. Thus we may write 


[Pa bf an =F ari). (2) 


Differentiating with respect to ¢, we find a differential equation, which, 
solved for r, gives 

r=2V/Dt. (3) 
Hence from (1) and (3) 


ome a (1 -rum): (4) 


If we put c(x,t) = }co, we find the relation between x and ¢ when the con- 
centration is at its half value. This relation is of practical use in the deter- 
mination of D. Solving for this particular value of x, namely, x*, as a 
function of ¢ we see that 


x*=/Di. (5) 
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This may be compared with the relation given by the exact solution (cf. 
Carslaw, 1945) in which «* = 0.9554/Dj. Similarly, for small values of 
x/+/4Dt, (co — c)/cy in the exact case has the value 0.56«/+/Dt, whereas 
the approximation method gives 0.5x/+/Dt. 

In some problems we might want to know the total amount of material 
which has flowed across « = 0. The approximate value given by (2) and 
(3) is seen to be co\/Di, whereas the exact value is 1.13¢o\/Di. 

To illustrate the effect of a different approximation for the relation 
between c and «x, we consider an exponential relation in which the dis- 
tance 7;(¢) is analogous to r(#) above: 

c (x,t) = cen), (6) 


Since 7;(0) = 0, therefore c(x,0) = 0 for x > 0; and unless 7(t) = 0, 
c(0,t) -— Cos 
Following the argument used to obtain expression (2), we find 


if Fat= Df 2 = Cor (t) , (7) 
0 0 zi 


the last term being the integral of the concentration from x = Otox = o. 
Differentiating and solving for 7;(¢) we find 
n= v2D, (8) 
so that ¢ is given by = 
c(x, t) = Gye */¥2D! (9) 


The relation x*(¢) for c(x*,t) = co/2 is then 
a*(t) = 0.980 VDt (10) 


as compared with 0.955+/ Dt. Similarly, for small values of «/+/2Dt, we find 
Co— C(%, bt) _ 0.707 x 
Co = V Di ¢ 
while the integral of the flow is co\/2Dt. Thus this approximation is not 
as good for small x/+/2Dt as the linear approximation, but is better for 
larger values. For many biological applications there would be little 
reason. to choose between these two approximations. 


Case 2. Same as Case 1 but for a finite region. 

Consider a finite region of length L, containing no sources, in which the 
concentration is maintained at a constant value at « = 0, starting at 
time ¢ = 0. Initially the concentration is zero in the region. For a time 
after ¢ = 0 we may use the same expression as in Case J. However, when r 
becomes equal to the length of the region, at some time i*, the approxi- 
mation will no longer hold. Thus for values of ¢ which are less than the 
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value {* = L2/4D, obtained by equating the right side of (3) to L, ex- 
pression (1) with r given by equation (3) is the approximate solution. 

If at x = L the concentration is maintained equal to zero, then, for 
t > t*, ris equal to L and hence 


¢(%, i) = «(1-F); 


that is, the steady state is reached at the time ¢=¢* = L?/4D. Atx = L/2 
and when ¢ = ?*, the exact value of c is 0.446cp instead of 0.5¢o. 

If at x = L there is a reflecting barrier instead of a sink, then the con- 
centration at x = L will not remain zero but increase. Suppose that the 
concentration distribution remains linear. The average concentration ¢ in 
the region is the value of c at x = L/2. The flow into the region is given 
by D times the concentration gradient 2(¢) — €) /L. The integral of this 
flow from ¢ = ¢* to t is equal to the change in the total material in the re- 
gion in the time interval ¢ — #*, i.e., GL — col /2. Hence ¢é may be found 
as a function of time, and since the relation is taken to be linear with dis- 
tance in this method, c(x,t) is then found to be 


(ey = «| 1-4 exp (5-77)]- Gib) 


If the region is divided into two compartments instead of one, a some- 
what better approximation can be obtained. 


Case 3. One-dimensional semi-infinite or finite region with a source at 
x = 0. Initially c = 0 everywhere. 

Let there be a constant flow Fy into the region from the left at x = 0 
beginning at time ¢ = 0 (Fig. 2). We suppose that nothing diffuses to the 
left of « = 0. Let the concentration be zero in the region at ¢ = 0. Let 
c(x,t) again be given by the linear approximation (1). Then we again have 
equation (2) but with co a function of time in this case since the concentra- 
tion is not maintained at a constant value. Differentiating with respect 
to time and eliminating either cy and ¢p or 7 and r from the resulting 
equations, we find r = 1/2Di and cy = rFo/D. We may suppose that the 
region is infinite or that it is finite with the concentration maintained 
equal to zero at « = L. The former is equivalent to the latter if L > o. 
Introducing the above values of r and cy into equation (1), we obtain 


Qt Fox <l 
Dita = Ee 


MO-De op 


D° 


Cc ts, t) =F) 
(12) 


Gey five 


bo 
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For small values of x/./4Di and for L > o, the exact solution is 


0.8 0F, eee = 


For x*(t) given by 
ES ae 


BOs ty.> 2 


and for the case in which L— , the approximation method gives 


x*(t) = »/Dt/2, whereas the exact solution gives «*(t) = 0.986./Di/2. 
For finite L, the exact method presents a more difficult problem, whereas 


Co (t) 


C(x, t) ——> 


n(t) 
x -> 


FIGURE 2 


the approximate solution can again be solved by inspection of equation 
(12) and found to be x* = »/D#/2, t < t*; 2* = L/2,¢ > i*. 


Case 4. One-dimensional finite region with diffusion coefficient D; and in 
which c(x,0) = co; elsewhere D = D,, c(x,0) = 0. No sources. 

Consider a region (Fig. 3) of width 2R in which the diffusion coefficient 
is D; and in which there is initially a uniform concentration co. Let the 
midpoint of this region be the origin from which any distance p is meas- 
ured. On either side of this region let the diffusion coefficient be D, and 
let the concentration be zero initially. After a short time the concentra- 
tion just beyond R will rise whereas at distances just less than R the con- 
centration will fall. Near the center of the middle region the concentration 
will not change much until some time 7. Hence the concentration may be 
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approximated, as shown in Figure 3, by the continuous function defined 
for t < t by: 
c=¢ for. p= /x-=R- Fs 


=f p li hoa ee 


S (13) 
e=fa(1-£-*) for ~ “RS Sie 


and c = 0 elsewhere. The quantities 7, s, and f remain to be determined 
as functions of t. As long as s < R(t < 2) we have the following relations. 
The material lost from the region p = |x| < R must equal that gained 


e) R 


5 ee 


FIGURE 3 


in the region p > R. Similarly, the flow toward the boundary p = R must 
equal that away from the boundary. Thus for ¢ < 7 we may write 


(=f) sey=rfeo=2f (2272) ar, (14) 
and 
Sif Doel (15) 
from which r/s = »/D,/D,, and, for t < i, 
r=2/D,t, (17) 


S=2VDzit. (18) 
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wn 


Since s = R for ¢ = 1, we find from (18) that 7 = R?/4D,. Substitution 


of f, r, and s from equations (16), (17), and (18) into (13) gives the de- 
sired solution for c(p,t) for t < i. 


After the time ?/, the concentration at p = 0 can no longer be set ap- 
proximately equal to co. We shall denote the concentration at the origin 


by c*. Then for ¢ > 7, c(p,t) can be approximated by the continuous func- 
tion (Fig. 4) defined for ¢ > i by 


eg | HER 


whist (Rte) 


r 


R<pSR+r, c=0 R+r<p, (19) 


FIGURE 4 


where c* is a function of ¢. Note that c* is equal to cy for ¢ < 7. We may 
then obtain, as in the above cases, the following expressions: 


‘Df c* 
IR — fRo* Rot =forr=2 fat, (20) 


Di(i=f)c* _Difo* 


? ; (21) 


Differentiating (20) and eliminating c* and f from (20) and (21), we 
obtain an expression for ¢ as a function of r: 


12) De i 2D rR) dre 
Dea RVD,/D; Dy? +2D:Rr+RD, 


(22) 


2D,(4¢—t*) =r? — 
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For a given r > R\/D,/D; = r(#), t can be determined by (22) in which 

the integral can be written in closed form. The value of r determines 

c(x,t) since, from (20) and (21), for ¢ > #, 

its Dee $i 

cs D; i ot DR : 

+= 2R (Dir + DR) (24) 
D;r+2D;Rr+DR 


(23) 


f 


If D; >> D., Di << D., or D; = D., the above expression can be 
simplified considerably. For the case in which D; = D., the rate of flow 
out of the inner region (within + R) is given by the approximation meth- 
od as coD/2/Di for t < i, whereas the exact solution gives coD/ /nDt 
+ ....The time at which } of the original amount has flowed out of 
the region is ? = R?/4D as compared with 0.82R?/4D. 

It may be noted from equations (23) and (24) that c* and f are given 
in terms of r. Therefore, if the sequence of spatial distributions of con- 
centrations is desired without requiring a knowledge of the corresponding 
times, these concentration distributions may be readily calculated with r 
as a parameter since f and c* are given in terms of r by (23) and (24). 


Case 5. Two-dimensional flow with cylindrical symmetry. For p < R, 
c(p,0) = co. For p >R, c(p,0) = 0; D= o for p< R; D=D- for 
pea oR: 

Consider a two-dimensional region in which the diffusion coefficient is 
infinite within a circle of radius R but outside this circle it has a value D.. 
Initially the concentration is equal to cp within the circle of radius R but 
is zero outside this circle. After a short time the concentration ¢* within 
the circle will decrease, though it will be uniform throughout the region. 
Outside the region the concentration will fall toward zero, being approxi- 
mately zero at some distance 7, which increases with time. 

Let p be the distance from the center and let c(p,t) be the surface con- 
centration at p at the time #. Since D = © within the circle R, we have 
c(p,t) = c*(t) for p < R, with c*(0) = co. The flow out of this region is 
equal to the loss from the region and is also equal to that which is gained 
by the outside region. Thus if 


(p,1) = c#[ i Seo) for’ (RX pe Hye ORGS) 


Tf 
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r being a function of ¢ alone, we find 


t t IDPs ° Rt+r 
Jf Paton [ROE dt= aR (cy) c%) = 22 f e(p, tinde. (263 


From this we find 


- R PR yt a 
and 


apt 2r—1.5R log (1 +3 +5 —_ 
(28) 
4 VSET _ 3K tan (a at oye) 


Thus c(p,¢) is given parametrically since c* and ¢ can be calculated for 
any value of r from (27) and (28). For the value of ¢ thus obtained, ¢ can 
be calculated for any p in the range for which c is nonzero. If only a se- 
quence of spatial distributions of concentration is required, we may calcu- 
late the concentration from (25) as an algebraic function of p with r asa 
parameter by eliminating c* from (25) by means of equation (27). To ob- 
tain a closed expression for c(p,t) we may extrapolate the solution for 
small values of r to some other larger value. Note that at r = \/3R, the 
first and last terms in the denominator of (27) are equal. Therefore we 
will take that larger value to be \/3R, and then extrapolate the solution 
for large values of r back to \/3R and force continuity on c in the boun- 
dary conditions. We then obtain 


col 1 fe | pysk? 


NE) es TRADI o AD: 
oe ae ARE 
(p —R) 2 
fai 
“| eee 3R2 
c(p,t) = eooe D2?’ 1? Tp 
2.30-+——,—_ + 3,8 


c(p,t) being equal to c(R,t) for p < R, and c(p,t) being zero for p > R 
4+ 2./Di or p > 1.98R + Dt/R’, respectively, for the two regions. 
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In asimilar manner we arrive at an approximate expression for c(p,t) in 
the three-dimensional case with central symmetry: 


4R® col 1 _25 | a 
¢ (a!) = Geet TOR? / Di + 16R Di + 8 (Di) ‘S paD' 
(p —R) i 
AR col 1 — eR y DHE | 2 
c(p, 1) = 687De DB '- FAD: 


14.27R?+ 16.41R Di+ 


R Re 


In this case the function is divided at the point at which the constant 
term and the cubic term in the denominator of the expression for c* are 


equal. 


Case 6. One-, two-, or three-dimensional flow with central symmetry. 
D = constant. For p > ro, c(p,0) = 0, for p < roc(p,0)co(1 — p/ro). Sub- 
stance consumed at a rate, ac, proportional to concentration. 

Consider a region in which initially the concentration is very nearly 
zero everywhere except near the center, where it is greatest. We shall sup- 
pose that the initial distribution can be approximated satisfactorily by a 
linear function of distance p from the center, the concentration being ¢) at 
p = 0, and approximately zero beyond some distance 7». Let the sub- 
stance be consumed at a rate proportional to its own concentration, a be- 
ing the proportionality constant. For the case of production, a is negative. 
If é is the average concentration at time zero, then é) = ¢o/(n + 1), n 
being the dimensionality of the region. After a time ¢, the average con- 
centration é(¢) will have decreased while the concentration will be zero at 
some distance r(t) instead of at 7. Since the rate of loss of the substance is 
proportional to the total amount present, we obtain 

Pon a 
He) oe —atr, (31) 

In this case the relation between 7 and / can be found as follows. Con- 
sider the boundary p = r/2 at any time ¢ and consider, for simplicity, only 
the process of diffusion, that is, a = 0, although the result is unchanged 
if a # 0. The flow across this boundary in time At is equal to c(0,f)DAt/r 
times the appropriate area for a given m. But this is the difference between 
the amount to the right of r/2 at ¢ + At and the amount at #, that is, it 
is the difference between the space integral of 


c(0, £+At) [1 ~ al 
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taken from p = r(t)/2 to r(t + AZ), and the space integral of 


(0, 0) [1-8] 
2 
taken from r(t)/2 to r(t). The expression resulting from equating the dif- 
fusion flow to the difference of the above-mentioned space integrals, to- 
gether with the condition cr = constant, leads to the relation dr/dt 
= 4D/r which is independent of 7. For example, when x = 1 the resulting 
expression becomes 


37 (1) = 
ees = OO) 


ef bo gr \ Al 
=(r Gases) ay : 


But from the conservation of mass of the substance we have ér = Goro, 
so that rd¢/dt + cdr/dt = 0. Multiplying this by } and subtracting from 
the above equation we find dr/dt = 4D/r. This expression is found also 
if 2 is two or three, and if a ~ 0. Upon integrating we find 


2¢(t+ At) DAt _ 
r 


e(t-+at) v(t) [1+ 


i= V8 Di+ r?. S52) 


This expression satisfies the requirement that r “ roatt = 0. From equa- 
tions (31), (32), and (1), and noting that ¢(¢ = 0) = co/(m + 1), we ob- 
tain for c(p,t) the expression 


—at somes Sa BS ep 
ar (: V8 Di+ ze) 


n/2 i 
(+52) 
"9 


the value of c being zero wherever the expression leads to negative values. 

Expression (33) may be compared with the corresponding exact solu- 
tion for the case in which 7) approaches zero but in which coro remains 
finite. If OQ; = coro and Qn = mroco/3 for n = 2 or 3, then the exact solu- 
tion becomes (cf. Carslaw, Joc. cit.) 


eat ep? /4Dt 
c(p; #) = Ci 


c(p, 1) = (33) 


If ro is permitted to approach zero in expression (33), the resulting ex- 

pression can be seen to be a fairly satisfactory approximation to (34). 
Consider next a finite region with a reflecting barrier at p = L. Then 

equation (33) would be approximately correct until the time i* when 
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p = L. From (32) this time is given by * = (L? — 75) /8D. If we let ro ap- 
proach zero then for n = 1 and t > i* we find (Q; = 2éoro = Coro): 


c(p, ) =2t[1+(1—52) exp -)] ent, (35) 


In order to obtain (35) the region was divided into two compartments. 
Otherwise the method is the same as that used to obtain (11). The corre- 
sponding exact solution for this case with a reflecting barrier can be writ- 
ten as 


6(p; t) aa» { e—@nL—p)*/ADty eons . (36) 


n=1 


A comparison of this expression with the approximate solution from 
(32) and (34) again shows that the results of the approximation method 
are quite satisfactory. For example, the ratio of the values of c(L/2,i*) 
as given by the two methods is 0.985. 

If the initial distribution can be approximated by a rectangular dis- 
tribution instead of that used in the preceding paragraphs, then the meth- 
od of Case 4 can be followed to obtain c(p,t). If 4; is the time at which 
the distribution by this method becomes triangular (which happens for 
t, = i* of Case 4) then, for ¢ > 4, we again find expression (33) if ¢ is re- 
placed by ¢ — t:, ro is replaced by r(t), and if m = 1. 


Case 7. Same as Case 6, except rate of consumption, Bc, proportional to 
the square of the concentration. 


In this case the approximation (32) will be less satisfactory than in the 
previous case. However, even such a rough approximation can be ex- 
pected to show the effect of the various parameters on the temporal and 
spatial distribution of the concentration. 

Consider a situation which is the same as that of Case 6, but with the 
consumption proportional to the square of the concentration of the sub- 
stance, the proportionality coefficient being denoted by 8. If we now fol- 
low the method used to obtain (32), we find that the quantity 8 does not 
drop out. If 6 is small enough the effect of neglecting it will not be serious. 
Therefore we can use (32). This approximation, even for moderately 
large 8, can be expected to show the effect of the various parameters on 
the temporal and spatial distribution of the concentration. 

The rate of decrease in the total amount of material can be found by 


et 
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integrating Bc? over the whole space. We thus obtain, after canceling out 
numerical factors, 


ad (ér*) ae # 

Pee A oa yh das Me (37) 
dl B cr ’ B n(+2) . \o 4) 
Eliminating r from (32) and solving for ¢ as a function of f and then using 
equation (1) with « = p, we find the following approximate expression 
for the concentration as it depends upon distance and time: 


co (1— Seeqe GE =) 
eae (38) 


; B ae pen 
—n/2 
Tep@+1) SS 


The quantity in the denominator involves the logarithm for m = 2; other- 
wise it involves the same radical as appears in the numerator. 

If in expression (38) 2 is equal to 2 and the factor co/(1 — 8Dir-’) is 
moved to the left-hand side, the resulting expression may be compared 
with equation (37) of Landau (loc. cit.), which is an approximate solution 
of the same problem by a more elaborate method. The two expressions 
will be found to be the same if 6? = 79/2 and a = 36/2. Actually B should 
be equal to a, while r) should be equal to 1.8 6. Thus expression (38) is 
fairly satisfactory. 

If both monomolecular and bimolecular reactions occur, the denomina- 
tor becomes somewhat more complex and the factor e~** occurs in the 
numerator. 

The author is indebted to Dr. George Karreman and to Mr. C.S. Patlak 
for reading and discussing the manuscript. 
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The theory of imitative behavior, developed previously, is applied to the case of two social 
groups which are separated spatially. If the information of each group as to the behavior of the 
other is complete, the case reduces to that of a single group. When any information is lacking 
at all, the two groups are independent. If we have two mutually exclusive behaviors 4 and B, 
all four combinations AA, AB, BA, and BB are possible. If the mutual information gradually 
increases from zero, then for a certain value of it, the group which is more informed about the 
behavior of the other will change to that behavior if it did not already exhibit it. If for constant 
information the size of the group increases, then, above a certain threshold value, the larger 
group imposes its behavior on the smaller. 


In previous papers (Rashevsky, 1949a,b; 1950) and in our book 
(Rashevsky, 1951a, hereafter referred to as SB) we have studied various 
aspects of imitative behavior in large social groups. With one exception 
(Rashevsky, 1950) all the above studies applied to the case of a popula- 
tion distributed uniformly in space. The one exception introduced the 
possibility of a continuously varying population density. In reality, how- 
ever, we more frequently have almost discontinuously varying population 
densities. Human beings concentrate in more or less large groups such as 
cities, towns, or villages, where the population density is relatively high. 
In between these cities, towns, or villages the population density is very 
much smaller, and in many cases may be taken as zero with a sufficient 
degree of approximation. 

In the present paper we shall inquire how such a quasidiscontinuous 
distribution affects various aspects of imitative behavior, which may de- 
termine the adoption by the whole population of a uniform behavior pat- 
tern. Such a study may be of particular interest for the understanding of 
the earlier phases of human history, when the quasidiscontinuity of the 
population density was particularly pronounced. The city states of an- 
tiquity and their interaction may present possible applications of the 


study. 
63 


64 N. RASHEVSKY 


Consider first two social groups, or communities, which are spatially 
separated. Let the population of the first one be Nj, that of the other V2. 
We shall also again restrict ourselves to the case of only two behaviors, 
A and B, which are mutually exclusive. The notations used will be the 
same as in SB; the quantities referring to the two groups will be denoted 
by corresponding subscripts. We shall assume that the average preferences 
gi and ¢2 for A or B are zero for both groups: 


é1=¢2 = 0. (1) 


Positive @ shall be taken to correspond to a preference for A, negative 
for B. Here we use the notations A and B, instead of R; and R2 as in SB, 
because in our present notation different indices refer to different groups. 

If the two groups are completely isolated, so that individuals of one 
group do not know anything about the behavior of individuals of the 
other, then the fundamental equations of imitative behavior for each 
group are independent. We have 


7 

Vt = Ay (X1— Vi) = as (2) 
and 

d 

a= 2(X2— V2) — dee. (3) 


The rate of change of y, is determined only by the number of individuals 
X, and VY, in the first group who exhibit behaviors A and B correspond- 
ingly. A similar thing holds for dy2/dt. 

If, on the other hand, due to an exchange of information between the 
two groups, the individuals of each group are fully informed of the be- 
havior of the individuals of the other group, then we have 


d 

oe 1 PO a ee ee (4) 
and 

d 

OS = a CR ere eeie (5) 


In general, when the information is incomplete the effect of X, — VY. 
on the rate of change dy,/dt will be less than the effect of (X, — VY). 
Vice versa, the effect of X; — VY; on dy2/dt will be less than that of O& 
— Y2). Denoting by « and €: two quantities such that 


0OSaS1, 0Sesi; (6) 


— 
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we have 
d 
Hs = 1[ (Xi — Vi) +a (X2— Ve)] — ayn, (7) 
dy - i 
ap = As (a Yo) te(X1— Vy) ] — ar. (8) 
Putting 
X1— Yi = NF i (y1), (9) 
X2— Yo= NoF 2 (pe), (10) 


we obtain the nonlinear system: 


eae Ai Nii (¥1) +a AiNMoPs (Yo) — ani, C11) 
oY? 2A, WGP I 
77 tin te 2 NoFs (W2) +6 A2 Ni, (Y1) — dhe. (12) 


For & = € = 0 the system (11)—(12) reduces to two independent equa- 
tions. An analytical study of the system for any « and e; appears very 
difficult. We shall study the properties of the system geometrically. 

The requirement 


dy, 
= 50 (13) 


leads, because of (11), to the requirement 
€:A1NoF2 (2) 2 ai — Ai NF (1). Oey 


We shall denote the left side of (14) by L(2), the right side by R(y1). The 
function R(y2) is represented by the dotted line in Figure 1. Let R,,1 be 
the absolute value of the maxima and minima of R(y). 
The function L(y2) is represented by a sigmoid curve with asymptotes 
+1, with L(0) = 0, and 
dL (2) 
dps 


Therefore the inverse function ¥2(L) is represented by the curve of Figure 
2= Kor i 6 AiNo, v2 = +o, and fon —«A,N2, Yo =>—o, 
Consider the case of the equality sign in expression (14). If 


>0. (15) 


itn oe (16) 


m1 
A,N,’ 
then 2 becomes positively infinite for such values of ¥; for which R(y) 
= «A,N2 < Rm. For such values of y for which R(¥1) = —4AiN2, the 
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value of 2. becomes negatively infinite. Because Y.2=0 for L=0 
W2 is equal to zero for such values of y for which Ri) = 0. Therelote 
the equality sign in expression (14) leads to a curve in the Wi, W2 plane 
which consists of three separate branches, as shown by the full lines s 
Figure 1. Because of inequality (15) the inequality sign of (14) holds at 
all points “above”’ the three curves. This means that for all points in 
region I of Figure 1 we have dy,/dt > 0, while for all points in region I] 
of Figure 1 we have dy./dt < 0. On the three branches of the curve 
dy2/dt = 0. 
Because of (12) the requirement 
dy» 
dt 


= 0 Ge) 
leads to 

@ Ao N1F, (y) > Ao. — Ao Nok» (Wo). (18) 
Denoting by R,.2 the absolute value of the maxima or minima of a2» 
— A,N2F2(y2), we find, in a similar way as before, for 


Rio 
A2N, 


e < (19) 
the graph of expression (18) for the case of the equality sign in the yx, Yr 
plane. The curves are similar to those shown in Figure 1. They now run 
horizontally, however, since y and y2 are interchanged. Figure 3 shows 
in the yi, ¥2 plane both the lines of dy,/dt = 0 (full lines) and dy2/di = 0 
(broken lines). 

In the same manner we used in determining the sign of yi = dy /dt 
in Figure 1 we determine the sign of y2 in Figure 3. 

The nine points of intersection in Figure 3, O, Oi, O2, O3, O4, Os, Or; 
O;, and Os, represent configurations of equilibrium. From the signs of 
yn and yz in the different regions of the yi, ¥2 plane we can find the direc- 
tion of motion of a representative point in the neighborhood of any of the 
equilibria. These directions are indicated by arrows. We see that only the 
four equilibrium configurations represented by the points Oy, Os, Os, and 
O; are stable. Point O; corresponds to the case in which the first group 
shows behavior B predominantly, the second A. Point O; corresponds to 
the case in which both groups show behavior A predominantly; point O; 
corresponds to the case where the first group shows behavior A predomi- 
nantly, the second behavior B. Finally, point O; corresponds to both 
groups showing behavior B predominantly. Thus all four combinations, 
AA, BB, AB, and BA, are possible as stable configurations. We notice, 


a 
$,<0; Y,<0 


<0 
> b<o 
ay g,>0 
Sy are 
ed es % <0; f>0 
ee) 


FIGURE 3 


<0; p<o 


FIcure 4 
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however, that for 4A and BB (points O; and O;) the ordinates y, and yy» 
are greater than for AB or BA (points O; and O;). This means that if 
both groups show the same behavior, the majority of those who accept 
the preponderant behavior is greater than if each group shows a different 
preponderant behavior. 

If « and 2 decrease, the lines ¥; = 0 and yo = 0 approach straight 
lines, as is readily seen from Figure 2. In the limit e, = & = 0 we have 
the situation shown in Figure 4. All four stable configurations are now 
symmetrically situated. That is, whether we have AA or BB, or AB or 
BA, the majority exhibiting A (or B) in the first group is always the same, 
regardless of what the other group exhibits. The same holds mutatis mu- 
tandts for the second group. This is what the case should be if the two 
groups are completely independent. It is also important to notice that for 
€, = €2 = 0 the majorities in the cases AA and BB are smaller than in 
the previous case. This is seen from the fact that the coordinates of O; 
and O; are larger in the case represented in Figure 3 than in the case 
shown in Figure 4. 

What happens, however, if € and €2 increase? As « increases the points 
O; and O2 move closer together, as well as do the points O; and O,. For a 
certain critical value «* of 4, the points O, and O, will fuse; at the same 
time O; and O, will fuse. The resulting points, as is readily seen, will be 
unstable. Moreover it is seen that the instability will be such that if the 
system was in configuration O, it will move to O;; if in configuration O; it 
will move to O;. Hence we find the following result. 

If for a constant €, the value of e increases and reaches or exceeds a 
critical value e,*, which is a function of €2, then the combinations AB and 
BA become unstable. If the system was in configuration BA it will pass 
into AA; if it was in AB it will pass into BB. The same holds when « is 
kept constant while e, increases. However, now the configuration BA will 
pass into BB, while AB will pass into AA. 

We thus see that for values of ¢ and €2 which exceed certain thresholds 
the two spatially distant groups will behave as one group. They will 
either both exhibit behavior A or behavior B. This ‘‘complete coupling”’ 
will occur when e, becomes sufficiently larger than Rmi/AiN2 and/or €2 
becomes sufficiently larger than Rn2/AoMi1. 

Since €, enters in the combination «Ai, and e enters in the combina- 
tion €2A2Nj, it follows that instead of increasing « we may increase Ns, 
or instead of increasing ¢2 we may increase N;. We have seen that an in- 
crease of ¢ results in the imposition of the behavior of the second group 
on the first. Hence an increase of N» will also result in an imposition of the 
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behavior of the second group on the first, while an increase of JN, will 
result in the second group adopting the behavior of the first. The larger 


of the two groups predominates. 
A very crude evaluation of the critical magnitude of ¢ and €, may be 


made if we approximate F\(y) and F2(y) by 
F(y) =1—e-. (20) 


Let us also put, for simplicity, Fi(y) = Fo(y) = FW), Ai = Ao = A, 
a, = dy = a, & = €& = e. Consider the behavioral combination AB. 


FiGurE 5 


When ¢ = 0 the behavioral configuration of the first group is deter- 
mined by the intersection of the curve AN,F(y) and the straight line 
ay, since we have 


ANF (vi) — ayi=0. (21) 
When 0 < e < 1, then we have for the first group, instead of (21), 
A NiF (f1) —€A NoF (Y2) — ay. = 0. (22) 


The transition from (21) to (22) is formally equivalent to a vertical 
shift of the curve AN,F(y,) downward by an amount 6 = eAN2F(y2). 
As ¢ increases the shift becomes greater until the situation shown in Fig- 
ure 5 is reached. Then the behavior of the first group changes suddenly 
from A to B. 

In the position shown in Figure 5 the equation of the curve is 


ANF (yi) = AN, (1 — e*h1) — 8. (23) 
At the point of tangency, whose abscissa is y*, we have 


AN (1— e-™) — 5 = ay* (24) 
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and also 
A Nik e** =a, (25) 


From (25) we have 


Sees | 
¥ ee 


Equations (24), (25), and (26) give 


ee at, (26) 


ES a a a 
eS we ig 
ey te Be ANe One 
or, putting a/kAN; =w <1: 
§6=ANi(1—w+w log w). (28) 
Since 
6=eA NoF (Wo) ) (29) 
therefore 
al See +w lo ) (30) 
F (v2) N2 ® @ g i). 


The value of F(W2) is of the order of magnitude of one. The value of w, 
and therefore of the expression in parentheses, is usually also of that order 
of magnitude unless we consider extremely small values of a. Hence, in 
order for the first group to change from behavior A to B when the second 
group exhibits B, e must be of the order of N;/Ne. By a similar argument 
we find that if the second group is to change from B to A, if the first ex- 
hibits A, e should be of the order of N2/N;. If N; > No and e€ increases 
from zero, it will first reach the value N2/N;, and the second group will 
follow the behavior of the first. If Nz > Ni, then as soon as e€ > Ni/N2, 
the first group will change from A to B. In other words, the group with 
the larger population eventually imposes its behavior pattern upon the 
other group, as we have already found from geometric considerations. Or, 
if both groups are of the same size, then the group which is better known 
to the other imposes its behavior. 
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It is demonstrated that, if the variations of viscosity throughout a cell are considered, swell- 
ing stresses may produce elongation and division. To do this it is necessary to generalize Betti’s 
theorem to cover systems containing viscosity gradients and such a generalization is presented. 

On the basis of two special assumptions it is shown that most of the results of the diffusion 
drag theory of cell division may be duplicated by the present theory. 


I. Introduction. It is well known that many substances dissolved in a 
gel may change the specific volume of the gel. If the concentration of sub- 
stances is not uniform, stresses will be produced. These swelling stresses 
will, in general, produce flow. N. Rashevsky (1948a) investigated such 
stresses to see if they could produce cell division. His conclusion was that 
swelling stresses could not produce elongation and hence could not lead to 
cell division. 

In this paper a slightly more complex model of a cell than that analyzed 
by Rashevsky will be considered. Whereas Rashevsky studied models of 
uniform viscosity, in our model regional variations of viscosity will be con- 
sidered. It will be found that this one additional feature will have an impor- 
tant bearing on the question of whether or not a cell can elongate under the 
influence of swelling stresses. It will be found that cells containing viscosity 
gradients can elongate under swelling stresses even though cells of uniform 
viscosity cannot. 

One result of this analysis of elongation by swelling stresses will be a 
qualitative explanation of the inhibition of cell division by applied hydro- 
static pressure (Marsland, 1942). ; 

II. The generalized Betti theorem. Betti’s theorem has played a promi- 
nent role in all theoretical studies of cell division (Rashevsky, 1948b). This 
theorem yields the average strain rate in a hydrodynamic system in terms 
of certain integrals involving the body and surface forces acting on the 
object under investigation. It therefore enables one to obtain a measure of 
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the elongation of a body without solving the extremely difficult differential 
equations of hydrodynamics. 

To date, Betti’s theorem has been derived and used only for bodies of 
constant viscosity. The derivation has proceeded by way of a reciprocity 
relation in elasticity theory and an analogy between elastic systems and 
incompressible fluids (Young, 1939). The present work demands a generali- 
zation of Betti’s theorem, one valid for systems containing viscosity gra- 
dients. As a step leading to this generalization a new derivation of the 
classical Betti theorem will be presented. This derivation is simple and will 
immediately indicate the means of generalization. It will, incidentally, re- 
move the restriction of incompressibility and demonstrate the form of the 
theorem for compressible fluids. 

Neglecting inertial forces, the equation of equilibrium, both in elasticity 
and hydrodynamics, is* 

003; 
Ox; 


ae et) (1) 


where g;; is the stress tensor, f; is the body force (per unit volume), and x; 
is the ith Cartesian coordinate. 
The boundary condition is 


X 5 = 055"; ,; (2) 


where X; is the surface force and v; is a unit vector normal to the surface. 
Let us denote by J the integral of the stress tensor over the volume of the 
body. Then 


T=JSffoiidvV =JSSSoindjdV ’ 
where 6;; is the Kroneker delta. 
Since 


petty 
Ox,’ 


ce] Ox; 
Pa lSiexs e¥ 


5x; 


* Geli : 
Tensor notation is used throughout this paper. The summation convention is assumed, 


ae two identical indices in any term imply that the term is summed over all values of the index. 
us 


Xi Gis =%1Qig + Xo gay + 393; and 
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and, integrating by parts, 


T= SfounesdS — ff fx; 52% av. 
» k 


Therefore, using equations (1) and (2), 


Si feadV = ffaxj,X.dS+ fffasadv. (3) 


Equation (3) is completely equivalent to Betti’s theorem as may be 
shown by applying it to the fundamental equation of Newtonian hydro- 
dynamics (Prandtl and Tietjens, 1934, p. 256 ff.) 


it : 1 
ot rte | Or Oi3 = 2p (ej —3 €k 0:3) ’ (4) 


u being the viscosity and e;; the strain rate. 
Integrating (4) over the volume of the cell and remembering that 


: d 
SSfendv ==, 
we obtain 
1 dV Ree a? Ta oe, 
SSfesdV =3 Pe Jet ae Ji eiX eX a 
1 
5 SS Senndisd V = 5 SSeehedisd 
A special case of (5) is 
1 
SffedV =3 +z SSS fai (whataahs) {AV 
(6) 


1 


which, apart from the first term on the right-hand side, is a familiar form in 
the literature of mathematical biophysics. 

It should be noted that Betti’s theorem gives us no information about the 
over-all dilatation of a system. The integrals in equations (5) and (6) ex- 
press only the distortional part of the over-all flow. oa 

The means of generalizing Betti’s theorem is now clear. If our aim is to 
obtain an expression for 


SS SewdV 


76 IRVIN ISENBERG 


then, examining equation (4), we must evaluate J, which is defined by 


aS ffeoudV. 


Proceeding as before 
1 


Je did panels 
=S If vn si dV 


= Sf = ou jaS — SIs x; Sat av 


Ox 
- 116000 2,(2)¢ ie 


and, therefore, using equations (1) and (2), 


en 1 it dfs 
Tees dV = SS > XimydX + SSS = fs V — SS fave ge (S) av. 


Therefore 


1dV 1 ial ok e7ece = 
SffegdV =z di bata lis Maisto ISS des 


1 a) 1 

—F SS Sours s-(4) av 
1 1 il 1 

— 6 fe ares SSS = fatadisdV 
1 07! 

Te ise uti g-(-) 6,;;dV. 


This is the required theorem. 

It should be noted that this theorem is much more difficult to apply than 
the classical Betti theorem. For an actual evaluation we must know not 
only the surface and body forces but the stresses as well. This will, in gen- 
eral, require the solution of the set of equations (1), (2), and (4). But, if one 
could solve these equations one would have no need of Betti’s theorem! 
However, certain qualitative conclusions may be reached by the use of 
equation (7). 

ITT. The equivalence theorem for swelling stresses. Let us consider a cell 
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in which there are swelling stresses but no body or surface forces. The flow 
in such a cell would be governed [see equations (1), (2), and (4)] by 


Ooyj; 


aa (8) 
o;v;=0, (9) 
1 : ibe, 
od Sars and si = 2a( is — 3 nd : (10) 
K 
enero Ore +y (c) . (ll) 


The meaning of equation (11) is as follows: o;:/3 is the hydrostatic ten- 
sion at a point in the cell; K is the compressibility, while e;; is the dilata- 
tion; (c) would be the dilatation of a volume element if it were removed 
from the cell and had a concentration c of whatever substance is producing 
the swelling. Thus equation (11) means that the dilatation at any point is 
the superposition of the free expansion and that due to the existing 
hydrostatic tension. 

Suppose we define o;; by 


— 1 
Tig = Fis — HOG: (12) 


Then equations (8) through (11) become 


doi _1 oY _g (13) 
Oe KK OF; : 

ania (14) 

Gs Mendes = 2m ees ods : (i) 

a (16) 


€kk <a y Ckk - 


Thus the flow in a system with swelling stresses is the same as the flow in 
a system with no swelling stresses but with applied forces 


1 oy 
(i eteg ae 
K 0x; (17) 
1 
X:=sW 
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It is worthy of note that the above equivalence theorem is the same as 
the theorem concerning strains in an elastic body subjected to thermal 
stress (Love, 1944). It is easy to see why this should be so. Equations (8), 
(9), and (11) are the same for hydrodynamics and elasticity, while equation 
(10) is invariant under the transformation o;;—> oi. 

The above equivalence theorem will not be used in the remainder of this 
paper. In spite of this, however, there is some justification for its inclusion 
here. It seems rather difficult to directly visualize the effects of swelling 
stresses on elongation. It is much easier to discuss the effects of applied 
forces on the cell. The above theorem demonstrates that this is possible and 
equation (17) gives the form of the forces which may be presumed to act on 
the cell. Furthermore, the theorem demonstrates that the so-called “con- 
tradictory”’ assumption of Rashevsky (1948a), whereby the stresses are 
calculated as if the cell were a solid while the flow is treated as if the cell 
were a viscous liquid, is not really necessary. 

IV. Cell division by swelling stresses. With no body or surface forces the 
generalized Betti theorem becomes 


1dV 1 0/1 : 
SffesdV = 3 oat oe era 
(18) 
te 5 SS fewe: ~ (=) av ius, 
a special case of which is 
io Vo 
feud Vg a a Ss ee ~(<)a V 
(19) 


1 0 1 
+5 SSSeur 5— (=) av. 


The above integrals will not, in general, vanish. If the viscosity is the 
same throughout the cell, however, 


Sffedv =4 Y, 
and there will be no elongation. This duplicates Rashevsky’s result (1948a). 
It can be seen that the only stresses important in elongation are those 
at regions where gradients of viscosity exist. The stresses at other places 
in the cell will give rise to streaming but will not contribute to the over-all 
elongation of the cell. Since this is the case one might suspect that the most 
important stresses will be those at regions where the gradients of viscosity 
are greatest. Such a region is the inner surface of the cortical plasmagel 
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where the viscosity changes abruptly. Indeed , during mitosis, while viscos- 
ity changes occur throughout the cell, the increase of the consistency of the 
periphery is quite marked (Seifriz, 1952), thus establishing a larger differ- 
ence in viscosity of the plasmagel and plasmasol than is present in a resting 
cell. 

If we do assume that only stresses at the inner surface of the cortex are 
significant, then equation (19) becomes 


1dV 
SSSfendvV =3—-> —-=) f fournds 


where jy is the viscosity of the outer region, y, is the viscosity of the inner, 
and S is the bounding surface. 

It is now possible to understand, on the present picture, the work on the 
inhibition of cell division by applied hydrostatic pressure (Marsland, 
1942). 

Protoplasmic gels belong to that class of gels which solate under pres- 
sure. Consequently yw, should decrease with pressure and 


(20) 


should become smaller. When this factor becomes sufficiently small the 
last two terms in equation (20) will be just balanced by surface tension 
forces, and cell division will be completely inhibited. 

Throughout this paper it has been assumed that a gel could be con- 
sidered to be a very viscous liquid. This is, of course, not true. Gels are 
solids which will, however, yield under sufficient shear stress. To treat the 
gel as a viscous liquid means that a model is proposed in which it is as- 
sumed (or hoped) that the “solid” properties of the gel play only a minor 
role in the phenomena discussed. 

If two further assumptions are made, equation (20) will yield almost all 
of the results of the diffusion drag theory (Rashevsky, 1948b). Since on 
one side of the surface S is the relatively fluid plasmasol it does not seem 
unreasonable that there should be no shear stress across ‘this surface. 


Thus, at the surface, 
(21) 


O55 = TEV VG 


and, neglecting the term 
av 
dt 


| 


80 IRVIN ISENBERG 


and recalling that ps > ba, 


; i} 1 
SfJend V = 7 Jf fonneds——— fi fomorme-+nee) dS. (22) 


Let us also assume that the interaction among the various parts of the 
plasmagel are such that the actual dilatation of an element of plasmagel 
near the inner surface is small. Then, using equation (11), 


om — 2 (0) (23) 
and 


(farre 5 Sf f0 co fnas—3 (rtetwes) dS. (24) 


According to E. Posnjak (1912) ¥(c) will typically have the form 
y(c) =k'c" (25) 
for a number of gels and dissolved substances. However, equation (25) has 
no theoretical foundation but is an empirical fit of data. Consequently, in 
any work in which the derivatives of y play no role the data may be pre- 
sumed fitted, over the range of c which is significant, by a straight line 


v (c) =ke-+constant , (26) 
so that 


iene = gt J foefnm —5 (axa H09ts) tas. (27) 


This has the same form as the equation governing cell division in the 
diffusion drag theory (Rashevsky, 1948b) if, in that theory, the cell mem- 
brane is assumed to have an infinitely large permeability. Consequently 
equation (27) will duplicate the results of the diffusion drag theory when- 
ever the permeability is assumed to be large in the latter theory. This 
covers most of the results of that theory. The fact that the integral in 
equation (27) is over the interface between the plasmagel and plasmasol 
while in the diffusion drag force theory the surface integral is over the 
surface of the eggs makes very little difference in comparing the two 
theories since the plasmagel is quite thin. 

However, it should be emphasized that though equation (21) is prob- 
ably reasonable and the use of a linear relation for y(c) can be justified, the 
assumption embodied in equation (23) may not be a good one. 
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It is shown that a slight modification of a model of excitatory phenomena in irritable tissues, 
which has been treated before, exhibits spontaneous oscillations. The frequency of these 
oscillations and the time-course of the potential across the model membrane have been de- 
termined, together with the dependence of some of their characteristics on some important 
parameters, particularly (Ca**). 


Rhythmic discharges in living organisms have often been compared 
with nonlinear or relaxation oscillations. A. Fessard (1931) stresses the 
analogy between neural rhythms and relaxation oscillations. Also A. 
Arvanataki (1939) has shown that decalcified nerves and myocardium may 
show self-excited oscillations. B. van der Pol and J. M. van der Mark 
(1928) have observed several similarities between the rhythmic heartbeat 
and relaxation oscillations as, e.g., heart block and demultiplication of 
frequency. The model experiments of A. Bethe (1943), O. Schafer (1943), 
and A. Bethe and H. Shafer (1947) have further emphasized these simi- 
larities. F. Brink e¢ al. (1946) have shown that the average frequency of 
impulses in the squid giant axon depends on (Ca++). A. L. Hodgkin (1948) 
has obtained repetitive discharges in the Carcinus axon on stimulation 
with direct current. 

It is the purpose of the present paper to account for spontaneous dis- 
charges from an extension of a model of excitatory phenomena which has 
been recently treated (Karreman, 1951, hereafter referred to as K) . Al- 
though several interesting features of excitatory phenomena were obtained 
from that model, it was not possible to show that the model could exhibit 
repetitive discharges even if oscillations did not appear to be unlikely. In 
this work we shall study a slight modification of the same model. 

Consider the model of a membrane, outlined in Figure 1, in which the 
average concentration of the K+-ion in one layer is x and y in the other. 
Assume that the permeability of the inner boundary J is determined by 
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considerations similar to those which led to equation (58) of K so that we 
may write 
7 hy = Bra; (1) 

similarly for the boundary /J between the two layers we shall write 
hi=Bruy - (2) 


As indicated in Figure 1, we further suppose that the inside concentration 
K, of K+-ions is high compared with that, K2, on the outside, as is usually 
the case in irritable tissues. We will denote the thickness of these layers 
by \ and uw and that of the external diffusion layer by v. As indicated in 
Figure 1 we will neglect in this paper, for reasons of mathematical sim- 
plicity, the influence of diffusion of the K*-ion on the inside. 


<——=— } >< & i ae LEVEL 


OF (K*) 


——— MEMBRANE 


FIGURE 1 


In the usual way (cf. K) we arrive at the following differential equa- 
tions for the determination of x and y as functions of £: 


dx h ss h 

FES } et Geos, (3) 
dy hy ie - 

ree “—y pe (4) 


Let us normalize the concentrations as follows: 


x=xK2, (5) 
y= 9Ko, (6) 
Ki= KiK. . (7) 


Introducing (1), (2), (5), (6), and (7) into (3) and (4) we obtain, after 
some simple rearrangements, 


d Es 
ap = ot pa? (Ki — 8) — 32 (2-5) }, (8) 


dy_»X 
Ging lF @=9) —r G1}, (9) 
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with 
on, (10) 
p= a, (11). 
es G12) 


For mathematical simplicity we will assume \ = yu. Furthermore we will 
define t = at and in the rest of the paper we will simplify the notation by 
writing x for *, y for 9, z for , and ¢ for r. We will consider only the case 
for which K; = 50. Equations (8) and (9) then assume the form: 


ot = px? (S0—x) —y2(x—y) , (13) 
Pye Ls 14 
ay ey) — 7 1) (14) 


—y=0 
A” AS INDICATED 
—--x=0 


p=1/30 


FIGURE 2 


The shapes of the curves described by equations dx/di = 0 and dy/ dt 
= (0 are shown in Figure 2. The influence of changes in p on the curves 
with equation dx/dt = 0 is much less pronounced than that of changes 
in r on the curves with equation dy/d¢ = 0. For this reason and because 
we are interested in investigating the consequences of changes in the 
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parameters which change the ’s, we will simplify the treatment still more 
by assuming that p is constant (= 1/30) so that the curve described by 
dx/dt = 0 remains unchanged. Note that (Ca++) and AV enter through 
the 6’s and that since p depends upon the ratio of the B ’s, it is likely to be 
relatively independent of the concentration of Ca**-ions and applied 
potential. However, the main conclusions obtained in this paper retain 
their validity if this simplification is no longer made. 

The dependence of the properties of the solution of the set of equations 
(13) and (14) on a change of the parameter 7 is shown in Figure 2. But 


<2 


Ficure 3. Illustration of an integral curve of the system of equations (13) and (14) starting 
at P(30, 6) for the particular values p = 34, r = 150, and K,/Kz = 50. The final part of this 
curve is a limit cycle. 


equations (12) and (62) in K give the dependence of r on changes in 
(Ca++) and AV, r increasing with (Ca++) and with AV at the anode. For 
sufficiently large values of r there is only one stable equilibrium point in 
the physiological range (as is clearly illustrated in Figure 6). For a value 
of r within a certain range there are three equilibrium points, only one of 
which is stable. For still smaller values of r there is within the physiologi- 
cal range only one equilibrium point which corresponds to a stable oscilla- 
tion, as is illustrated in Figure 3. The trajectory starting at P approaches 
the closed limit cycle (Poincaré, 1892) and approaches it near Q; the part 
PQ corresponds to the transient response. Any trajectory such as this one 
corresponds to a spontaneous oscillation, for no matter where the starting 
point P is situated, except O, the point representing the values of the con- 
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centrations will move along it, going by itself. If it starts at O, only an 
infinitesimal disturbance is needed to start the oscillation. After the point 
has described the transient part which is dependent on the position of the 
starting point P, it will move along the limit cycle an indefinite number 
of times; therefore, the limit cycle corresponds to the steady state re- 
sponse consisting of repetitive discharges. By numerical integration along 
the limit cycle the period can be calculated in terms of the various param- 
eters. Assigning the following plausible values to the parameters: K, = 50 
X 10-° and K2 = 10-* in M/I,& = 2 X 10 sec.!, L = 10-4cm., and 
the thickness of the membrane \ + » = 10- cm., the period depends 
only on the value which is assigned to the parameter 8. The value used 
here for this parameter is 10°. It might be significant that a value of the 
same order of magnitude, namely, 3 X 10°, was used to obtain the agree- 
ment between the theoretically derived curve and the experimental points 
shown in Figure 20 of K. 

The time-course of the potential difference across the membrane can 
also be derived from the model above. As an illustration we will take for 
this difference 

Kige isael a 350 


i ae eer 


which is similar to the first term in (25) of K. Of course we might have 
taken instead 

F log : : 

Because x is close to 50 anyway, and moreover, if it is considered 
that the changes in x are relatively small (see K) the difference 
between the results obtained by the two formulae is rather small. An ac- 
tual comparison for r = 150 has been made in Figure 4. The time-course 
of the potential can be calculated for the closed cycle; it has been plotted 
for r = 60 in Figure 5. From this curve the amplitude @ and the slope S 
of the descending part of the time-course of the potential can be obtained. 
The frequency is easily derived from the value of the period. Calculating 
these quantities for several values of r we obtain them as functions of r. 
The results of these calculations are shown graphically in Figure 6. 

For the value 200 of r there is only one stable equilibrium point within 
the physiological range of concentrations (Fig. 7). Depending on the start- 
ing point, which is determined by the changes of y and x from their values 
in the stable equilibrium state, and therefore on the strength of the stimu- 
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Frcure 4. Comparison of the potential responses obtained from the mechanism for the two 
potential functions (RT/F) log (50/y) (solid) and (RT/F) log (x/y) (dotted), for which 
r = 60, p = 1/30. The scale for the time ¢(=7/2) is obtained by setting a = 10?. 
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Ficure 5. The same as for the solid curve in Figure 4 except that r = 60 


88 


» 
400+100 co 40 
200+ 50 ail 20 
A 8 
° Pee eee EE se eS ee ee EE ee es er ee ee ee i) 
oO 50 100 150 200 


nh (Cat *)—> 


FicurE 6. The frequency » (in sec.), the amplitude a (in mV.), and the slope s of the 
descending part of the voltage-time curve (in V. sec.) plotted as functions of r, which, in 
turn, is a function of Ca** concentration. 


FIGURE 7 
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lation, the trajectories have a different shape and size. For increasing 
strength of stimulation, and therefore increasing changes in y, the tra- 
jectories are ordered from J through VJ. It may be noted that there exists 
a threshold for the stimulus which has to be exceeded for an “active local 
response”’ to occur. The corresponding time-courses of the potential dif- 
ference across the membrane, obtained in a similar way as above, have 
been plotted in Figure 8. These also show the existence of a threshold for 
the stimulus; below a certain value of the strength of the stimulus the po- 
tential returns monotonically to its original value. Above a certain value 
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Ficure 8. The voltage, (RT/F) log (50/y), as function of time for several values of the 
strength of a suddenly-applied stimulus as calculated from equations (13) and (14), with 
a= 102, p= #0) r= 200. 


of the strength of the stimulus there occurs, after the negative artifact of 
the stimulus, an additional negative change in potential, which corre- 
sponds to some “active process.”’ After passing through a phase of “‘posi- 
tive after-potential” the potential returns to its original value. 

It might be useful to point out that this result would correspond to a 
local.response (cf. K) and therefore does not affect the all-or-none law; on 
the contrary it again (K) points out the possibility of the existence of a 
local threshold. 

The above results can account for a considerable number of phenomena 
of nerve behavior. They may also be applicable to the excitation of the 
heart. This is especially true in the case of spontaneous oscillations, so 
that this model may be useful as a starting point for theoretical investiga- 
tion of the heartbeat on the basis of physicochemical mechanisms. 
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Local regulation of blood flow as determined by capillary diameter and the number of open 
capillaries in a region is considered. The local changes in capillary diameter and in the number 
of open capillaries are assumed to be due to concentration changes of a diffusible, nonspecified 
metabolite. This metabolite is produced in the tissue and carried away by the blood stream. 

Using these assumptions and applying pertinent data on capillaries, deductions are made 
concerning: 

(a) the law of blood flow as a function of temperature and capillary radius for the hyperemia 
of high temperature, 

(b) high flow as it depends on metabolism during strenuous exercise of muscle, and 

(c) a first approximation to the time duration of occlusion hyperemia. 


Local regulation of blood flow by the capillaries has been observed by 
A. Krogh (1929) and M. van Dobben-Broekema and M. N. J. Dirken 
(1950). W. M. Bayliss (1902) and B. Folkow (1949) have presented evi- 
dence for a myogenic theory of regulation of capillary tone. This theory 
states that in accordance with known properties of smooth muscle, high 
capillary blood pressure is associated with capillary constriction, and low 
capillary pressure with capillary dilatation. The myogenic theory is 
utilized by Folkow (loc. cit.) to explain the transient effects of hyperemia. 
T. Lewis and R. Grant (1926) (also see reviews by Landis, 1934 and 
Pappenheimer, 1952) have given evidence for a local chemical action which 
causes dilatation or constriction. In this paper we will use the chemical 
action theory. We will apply recent data on capillary permeabilities to 
experiments on peripheral circulation in normal, hyperemic, and temper- 
ature-dependent flow conditions. 

If we assume a constant pressure drop across the capillaries, local capil- 
lary changes may change local blood flow by either or both of two mecha- 
nisms: 

(a) a change in the radii of the capillaries, 

(b) a change in the number of capillaries open to flow. 
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We assume the uniform production of a diffusible metabolite produced 
at the rate of g gm./cm.*/sec. throughout a homogeneous tissue. (The 
quantity g may be a function of the type of tissue,i .e., skin or muscle; the 
temperature, etc.) We further assume a uniform optimal concentration 
of this metabolite in the tissue, which we call Cr; (see Appendix I). The 
assumed mechanism of adjustment of blood flow is such that if the tissue 
concentration of metabolite rises above Cr the blood flow increases. This 
is accomplished by an increase in radii of the capillaries and/or an increase 
in the number of open capillaries. We assume the reaction to be caused by 
the raised tissue concentration of metabolite. Similarly, a decrease in con- 
centration below Cr is accompanied by a decrease of flow. 

Let us first consider (a) above: the mechanism of adjustment is a change 
in radii of the capillaries. This has been observed as the main mechanism 
of adjustment in the skin (see van Dobben-Broekema and Dirken, Joc. 
cit.). If R is the radius of the region supplying the capillary, / its length, 
and 7 the radius of the capillary itself, the metabolite produced in this 
region during a period of time At is r(R? — r*)lgAt. Since Cr is the steady 
state concentration of the metabolite in the region considered, the amount 
of metabolite diffusing across the capillary wall during the time Af? is 
2arlh(Cr — C.)At. Here h is the permeability of the capillary wall to the 
metabolite and C, is the concentration of metabolite in the capillary. Dur- 
ing the steady state the amount of metabolite produced in the tissue is 
equal to the amount carried away. Thus 


2arrlh (Cr—Coc) At= 7 (R?— r*) 1gAl. (1) 
We may simplify (1) by considering that r << R and assuming 
C. << Cr (Appendix IT). We then find: 
R 
r= Tics qd. (2) 
For the blood flow F, through a capillary we assume a law of the form 
Po =r, (3) 
where a will in general depend on the pressure difference across the capil- 
lary, the viscosity of the blood in the capillary, and the length of the 
capillary. 
For Poiseuille flow a = 4. 
Substituting (2) into (3) we obtain 


Re 
Fos 4 Fee g. (4) 
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In general q is dependent on the temperature. Assuming this dependence 
to be given by the Arrhenius equation, we have 


Gas Kk. B/Pr: (3) 
For the skin Qi9 = 1.8 (Krogh, Joc. cit., p. 274). Therefore we find 
B= 10,900 cal . 
Substituting (5) into (4), we find 
F=nF,= Kze7(8/RT)o | (6) 


with 2 equal to the number of capillaries in one cubic centimeter of tissue 
and F equal to the flow through one cubic centimeter of tissue. The re- 


LOGe VOL, FLOW 


a7°C 
Jer ef i a =o 


-00320 00325 -00330 -00335 0340 


! 
tr (AFTER FREEMAN) 


Ficure 1. Effect of changes in local temperature on rate of blood flow in normal hand. Or- 
dinate: Ine flow in cc./100cc. hand/min.; abscissae: 1/T (T in degrees Kelvin). 


sult, (6), might be compared with the data of N. E. Freeman (1935) for 
the skin (Fig. 1). The agreement in Figure 1 is obtained with the value 
2.8 of a, which is the right order of magnitude. The exponent a is smaller 
than would be expected on the basis of Poiseuille’s law. This might be 
due to the decreased viscosity of blood in small capillary tubes. 

We will now study the case in which the tissue responds to changing 
metabolic rate by a change in the number of capillaries draining the re- 
gion. This corresponds to mechanism (b) mentioned above. Suppose that 
there are capillaries of average length / per cubic centimeter of tissue. 
In equation (1) we considered the steady state condition in the region 
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drained by one capillary. As a modification of this equation, we now con- 
sider the production and drainage of metabolite in one cubic centimeter 
of tissue. Assuming, as before, r << Rand C, << Cy, we obtain 


2arinhC pAt = qAt (7) 


or z 
g=AhCr, (8) 


where A = 2zrlm is the total capillary wall area in one cubic centimeter 
of tissue. 

Since we assume that the pressure drop across the capillaries remains 
constant, the flow F through one cubic centimeter of tissue is directly 
proportional to the number of open capillaries. Thus, according to (7), 
the flow is directly proportional to the metabolic rate q. 

It has been observed (Krogh, loc. cit., p. 63) that many more capillaries 
of muscle are open during activity than are open during rest. In going from 
the resting state to the state of maximal activity the number of capil- 
laries open to blood flow may increase 40 to 50 times. According to J. F. 
Fulton (1949, p. 764) the volume flow during activity in muscle increases 
10 to 20 times. It is known (Fulton, Joc. cit., p. 811) that the oxygen utiliza- 
tion coefficient of active muscle is about 3 times that of resting muscle. 
The oxygen utilization coefficient is a measure of the amount of oxygen 
consumed per unit blood flow through the tissue. Therefore the increase 
in total oxygen consumption is expressed by the increased flow times the 
increased utilization coefficient. We thus see that the oxygen consumption 
of active muscle is about 30 to 60 times the resting value. In this case 
volume flow is not proportional to the number of open capillaries, pre- 
sumably because as many capillaries open the pressure drop across the 
capillaries decreases. However, we still find our expected proportionality 
between surface area of diffusion and rate of metabolism as expressed 
in (8). 

Observations on reactive hyperemia permit us to make an interesting 
application of our methods. Working with the reactive hyperemia of oc- 
clusion and strenuous muscular exercise, Lewis and Grant (Joc. cit.) con- 
cluded that this phenomenon is due to the local action of accumulated 
metabolic products. Equal initial blood inflow values were obtained from 
hyperemias due to (a) 10 minute occlusion and (b) 2 minute exercise of 
the forearm (Fig. 2, which is a plot of data of Lewis and Grant, loc. cit.). 
In the former case the hyperemia passed off in about 3 minutes and in the 
latter case in 10 to 14 minutes. 


In the steady state our tissue metabolite concentration Cy is maintained 
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by production of this metabolite at the rate q [see eq. (2)]. We assume 
that during the occlusion and exercise experiments the metabolite has 
been produced at a rate decreasing with time and has accumulated so 
that the tissue concentration of metabolite is greater than C r. On release 
of the occlusion this causes dilatation of the skin capillaries and the open- 
ing of muscle capillaries to flow. (The skin reddening observed with 10 
minute occlusion is caused by the dilatation mentioned above, whereas 
with the exercise hyperemia no reddening is noticed.) When the increased 
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Ficure 2, Comparison of hyperemia of 10 minute occlusion and 2 minute exercise. Inflow in 
cc./100cc. of tissue/min. (Lewis and Grant, Joc. cit.). 


flow has lowered the tissue concentration of metabolite to Cr the tissue 
returns to normal color with normal blood inflow rates. 

We will now apply our theory to obtain a first approximation to the 
above transient problem. If C(¢) is the concentration of metabolite in 
the region supplying the capillary at time ¢ after release of occlusion, then 
C(t) — Cr may be considered the excess concentration (above steady 
state value) at time ¢. The amount of metabolite flowing across the capil- 
lary wall during a small time interval Aé is then 


2arlh [C (t) —C,] At gms . 


The change in amount of substance in the region supplying the capillary 
during time At is r(R? — r)IA[C(#)] gms. The amount of substance pro- 
duced in the region during time At is 


a (R? — r?) lqgAt gms . 
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Thus by conservation considerations and simplifying by lettingr << R 
and C. << Cr (Appendix II), we find 

[2arlhC (t) — rR2lq] At = — tRIAC (2) . 
Making Af tend to zero, we obtain 


6G sATswollt 
2rhC—-R2q RB 


(9) 


Integrating (9) and applying the boundary conditions C(¢ = 0) = Co and 
C(t 0) = Cr = (R?/2rh)q, we find 
R? 


th f° (10) 


ie ee 
Et) = (c rar a) e @rh/R*)t 1 
The second term on the right side of (10) represents the steady state con- 
centration [see eq. (2)], and the first member the transient term. 

To evaluate C(#) from (10) we use average values of R and r. These 
are R = 100 X 10-‘cm. from Krogh (loc. cit.,p.62) andr = 5 X 10-*cm. 
from Fulton (loc. cit., p: 559). The quantity / is more difficult to derive 
and we use the results of a method suggested by the work of Pappenheimer 
et al. (loc. cit.). If we assume our metabolite molecule to be small and lipoid 
insoluble, then the permeability of the capillary wall will be determined 
by the fraction of wall available as “holes” for the diffusion of the mole- 
cule. We will call this fraction f. If D is the free diffusion coefficient of 
the molecule and 6 the thickness of the capillary wall, then 


D 
arek (11) 


We assume D to be equal to 10~° cm.’ sec.—', 6 = 10-4 cm. (Pappenheimer 
et al., loc. cit.), and f = 0.002 (Pappenheimer e¢ al., loc. cit.). Inserting 
these values into (11) we find that h = 2 X 10-* cm. sec.—. 

Using these values of R, 7, and # we will calculate C(#) from (10). The 
steady state value of C(¢) is given by 

2 
oy g=500q. 

We have assumed that the metabolite is drained from the region only 
by the blood. As a first approximation we assume a linear rate of accumu- 
lation during occlusion. Thus during occlusion it accumulates at the rate 
of q gms./cm.*/sec. Thus after an occlusion of 10 min., or 600 secs., the 
concentration of excess metabolite C. — Cr is 600g gms./cm?. This is of 
the same order of magnitude as the steady state value of C(t) which is 
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500g gms./cm*. Therefore, as a measure of the hyperemia duration we are 
justified in a choice of time equal to the time for the reduction of Cy — Cp 
to (1/e)(Co — Cr). We call this time t1,.. From (10) we find 


We evaluate ¢:,. from (12) using the above values of the parameters R, 
r, and kh. From this calculation we obtain the value #:;, = 8 minutes, 20 
seconds. This time is of the right order of magnitude considering the 3 
and 10 to 14 minute times observed as hyperemia duration (see Fig. 2). 


APPENDIX I 


An assumption implicit in the derivation of equation (1) and later 
equations is that the tissue concentration of metabolite Cr is uniform 
throughout the tissue. Using our values for tissue parameters we will now 
show that this assumption is valid. For the present we assume that Cr de- 
pends on the distance from the center of the capillary. We will call this 
radial distance x. Thus Cr = C7(z). 

The amount of metabolite carried across a single capillary wall during 
the time Af is (Cr(r) — C.)hAt/unit area. According to Fick’s law the 
amount carried to the wall during the time A¢ is 


dC 


See 


Al 


z>=r 


per unit area. From conservation considerations these are equal: 


dC ¢ 


[Cr (r) —C,] hAt = D ae 


Al. (13) 


rea 


In steady state conditions, the mass flow across any cylindrical surface 
coaxial with the capillary is equal to that across the capillary wall. Since 
this flow depends on the product of concentration gradient and area of 
surface, the greatest concentration gradient will exist at the wall of least 
area, or at the capillary wall itself. Thus 


dC 


C7 (R) —Cr(r) x ax 


Raye (14) 


To show that the concentration is approximately uniform throughout, it 
is sufficient to show that 


Cr (R) —Cr (7) 


iG: 
Cae ee ae 
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Substituting 
dC 
ix 2=F 
from (13) in (14) we find 
Cr (R) —Cr(r) SF ICr(r) —Cil (R= pin (16) 


Using the approximations that r << RandC. < < Cp we find from (16) 


Col(R) —Cr(r) = h Ra 


Cr(1r) Tra; oe 


Using our values for tissue parameters as in evaluating (11) we now find 


h R= 2G 1 Om ern sees 


D ~~ ~—«:10- cm.? sec. le 


And, from (17), 
Cr(R) —Cr( 


ut 5 
Crile 


which satisfies (15). 


APPENDIX II 


In our calculation we have assumed that C. may be neglected com- 
pared to Cr. We will now show that under conditions imposed by the in- 
creased blood flow of raised temperature and hyperemia this is justified. 
If the metabolite is adsorbed by some protein in the blood or undergoes a 
chemical reaction there, this condition will be satisfied. Otherwise we as- 
sume that it is removed before the blood re-enters the capillary (i.e., by 
the liver or kidney). Thus we assume that as the blood enters our idealized 
tissue the metabolite concentration C, is zero. 

The amount of metabolite diffusing into the blood in the time At is 
2rrh(Cr — C,)At/unit length of capillary. This quantity of substance is 
associated with a volume of blood of zr? cm.*/unit length of capillary. 
Thus the concentration of metabolite accumulated during time Af is 


2arh(Cr—C,) At 


rr 


gm./cm.’/ unit length of capillary . 
If we let C, equal the average metabolite concentration in the capillary, 
then 


2h (Cr—C,) 
“ ——— it 


2C,= 
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or 
2 C 
C,==- = (18) 
LCpry: 


where A+ is now the average time the blood spends in the capillary. Under 
increased flows of raised temperature or hyperemia we may assume that 
At is of the order of 1 second. Substituting this in (18) we find 


= 2 


Coma Cr. 


The quantity Cis thus seen to bea small part of C;. Assuming that @, ~ 0 
does not materially affect the results of this paper. 

The author is indebted to Drs. George Karreman and George Schmidt 
for their help in the preparation of this manuscript. 

This investigation was supported by a research grant H-627(C) from 
The National Heart Institute, of the National Institutes of Health, Pub- 
lic Health Service. 


LITERATURE 


Bayliss, W. M. 1902. ‘‘On the Local Reactions of the Arterial Wall to Changes of Internal 
Pressure.” Jour. Physiol., 28, 220-31. 
Dobben-Broekema, M. van and M. N. J. Dirken. 1950. ‘‘Reactions of the Vessels of the Rab- 
bit’s Ear in Response to Heating the Body.”’ Acta Physiol. Pharmacol. Neerl., 1, 562-83. 
Folkow, B. 1949. ‘‘Intravascular Pressure as a Factor Regulating the Tone of the Small Ves- 
sels.”” Acta Phys. Scandinav., 17, 289-310. 

Freeman, N. E. 1935. ‘“‘The Effect of Temperature on the Rate of Blood Flow in the Normal 
and in the Sympathectomized Hand.” Am. Jour. Physiol., 113, 384-98. 

Fulton, J. F. 1949. A Textbook of Physiology. Philadelphia and London: W. B. Saunders Com- 
pany. 

Krogh, A. 1929. The Anatomy and Physiology of Capillaries. New Haven: Yale University 
Press. 

Landis, E. M. 1934. ‘Capillary Pressure and Permeability.” Physiol. Rev., 14, 404-81. 

Lewis, T. and R. Grant. 1926. ‘‘Observations upon Reactive Hyperemia in Man.” Heart, 12, 
73-120. 

Pappenheimer, J. R. 1952. “Peripheral Circulation.” Ann. Rev. Physiol., 14, 259-82. 

., E. M. Renkin, and L. M. Borrero. 1951. ‘‘Filtration, Diffusion and Molecular Sieving 

through Peripheral Capillary Membranes.” Am. Jour. Physiol., 167, 13-46. 


. eu Mh, as 
U =a F oft ot arlene be. 
a ee Bee 


re an 


. 
® 


BULLETIN OF 
MATHEMATICAL BIOPHYSICS 
VOLUME 15, 1953 


SHORT NOTE 
NOTE ON A THEORY OF MEMBRANE PERMEABILITY 


RoBert A. WIJSMAN 


Division Or Mepicat Puysics, UNIVERSITY OF CALIFORNIA 
BERKELEY, CALIFORNIA 


In a series of papers (Bloch, 1944, 1946, 1947) I. Bloch has developed a 
model for the permeability of a membrane. The main feature of this 
model is that it associates with the membrane a potential function, de- 
scribing the forces between membrane molecules and solute molecule.The 
differential equation for the concentration of the solute is: 

‘pee 


+> =0 (1) 


ARTS 


in which c = concentration of the solute, V = potential function, 
D = diffusion constant, J = flow of solute across the membrane, k 
= Boltzmann’s constant, and T = absolute temperature. Here c, V, and 
D are functions of x, primes denote differentiation with respect to «, and 
the x-axis is perpendicular to the membrane. In Bloch’s theory T also is a 
function of x. In particular it was assumed that T(«) and V(x) are related 
approximately by: 


T(x) +350 (2) =T, = constant , 2a) 


which means that the total energy of a solute molecule does not vary while 
it passes through the membrane. 

The assumption, expressed in equation (2), seems very questionable in 
view of the collisions of the solute molecule with the membrane and sol- 
vent molecules, which tend to keep the average kinetic energy of the 
solute molecule constant. Equation (2) is an extreme case, and would be 
valid if the membrane thickness is small compared to the mean free path 
of the solute molecule in the membrane, so that collisions in the membrane 
can be neglected. In that case, however, equation (1) loses its validity. The 
other extreme case is that the membrane thickness is large compared to the 
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mean free path, and that c and V change slowly in a mean free path. In 
this case equation (1) can be applied, and in addition the temperature of 
the solute is approximately constant in x. The assumption T = constant 
seems therefore much better, and has the additional attraction of greater 
simplicity. 

After integration of equation (1) one obtains the following equation for 
the permeability 4 (defined as the ratio of J and the concentration differ- 
ence across the membrane): 


1 1 PSL). 
1 Se exp | RT | a, 33 


where the integration extends over the membrane region, i.e., the region 
where V(x) # 0. In the same way as in Bloch’s treatment (Bloch, 1944) 
equation (3) can be investigated after substituting expressions for D(«) 
and making assumptions on V(x). 
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BOOK REVIEW 


Dovuctas P. Murpuy. Heredity in Uterine Cancer. 1952. 128 pp. Cambridge, Massachusetts: 

Harvard University Press. $2.50. 

For many years medical research workers have attempted to study disease patterns in hu- 
man populations by sample survey techniques. Of the thousands of such studies which have 
been published, only a handful have provided useful results. The most common fault is an 
inadequacy of control data, although many studies have been biased by improper sampling 
methods. With the encouragement of The Commonwealth Fund which sponsors, edits, pro- 
duces, and distributes books and pamphlets as a philanthropic foundation, Dr. Murphy has 
produced a study of the occurrence of cancer among relatives of women suffering from cancer 
of the uterus. Aided in statistical matters of design and analysis by Dr. Lowell J. Reed and the 
late Charles P. Winsor, he has made a contribution to our understanding of hereditary factors 
in uterine cancer which does not suffer from the aforementioned deficiencies. For this reason 
the conclusions warrant acceptance and consideration by all investigators in cancer research. 

The book consists of four sections and a number of appendices. In Section I, a review of the 
literature reveals the inadequacy of previous information and justifies the need for the present 
study. The second section describes in great detail the methods used in the selection of pro- 
bands, the nature of the interrogation including follow-up studies, and the methods of handling 
the large amount of data collected. By far the largest portion of the book (about 50 per cent) 
is devoted to a tabulation of the data in Section III, according to many different criteria of 
classification. This section reads almost like a census report and exemplifies every experi- 
menter’s suppressed desire, so consistently thwarted by journal editors, to publish every scrap 
of data. It should be noted here that the readability of Section III would have been enhanced 
many fold if some of the more pertinent data had been presented in graphs. Section IV, only 9 
pages in length, contains a discussion of the results and the conclusions. Appendices A-E list 
a number of cooperating physicians, bureaus of vital statistics, hospitals, and other groups 
whose assistance is acknowledged. Appendices F and G reproduce the questionnaires and letter 
forms which were employed. 

The results of the study are based on the investigation of the female relatives of 201 cancer 
probands and 215 control probands. (A cancer proband is defined as an individual known to 
have suffered from uterine cancer and in whose family the occurrence of cancer is being 
studied. Control probands were selected at random without knowledge of whether they had 
previously suffered from cancer.) Information was obtained on 2809 relatives of cancer pro- 
bands and 3220 relatives of control probands. The conclusions were based on the family group 
as a unit which was defined as a proband, her mother and her sisters, her female first cousins, 
and their mothers. Beyond reasonable doubt it is established that the familial occurrence of 
uterine or non-uterine cancer predisposes a proband to uterine cancer, although the quanti- 
tative effect of heredity is small. On the other hand, no predisposition to cancer in other 
sites was found. By classifying the data in various ways, one can make many specific compari- 
sons but nothing much more definite can be said. In this connection it is noteworthy that the 
author does not attempt to make any conclusions beyond the scope of his study. 

For persons engaged in sample survey work of this nature, the book provides an excellent 
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example to follow. For the mathematical biologist, it points out a field worth investigating but . 
leaves much to be desired in the way of a foothold. Furthermore, the differences observed 
tween cancer and control family groups, though statistically significant at very low levels of it 
probability, are so small in magnitude that it seems doubtful whether a mathematical model 
sensitive enough to demonstrate such differences could ever be constructed. 


A. W. KIMBALL 
Oak Ridge National Laboratory 
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